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PART I 


STATISTICAL PREDICTION OF DYNAMIC 
DISTORTION OF INLET FLOW USING MINIMUM 
DYNAMIC MEASUREMENT — AN APPLICATION 
TO THE MELICK STATISTICAL METHOD 



SUMMARY 


A simplified explanation to the Melick method of inlet- 
flow dynamic distortion prediction by statistical means has. 
been included. A hypothetic vortex model is used as the 
basis of the mathematical formulations. The main variables 
of this model are identified by matching the theoretical 
total pressure rms ratio with the measured total pressure- 
rms ratio. Data comparisons using HiMAT inlet test data 
set indicate satisfactory prediction of the dynamic peak 
distortion for cases with boundary layer control device — - 
vortex generators. 

A method of the dynamic probe selection^ an essential 
part of this research, has been developed. Validity of the 
probe selection criteria has been demonstrated by comparing 
the reduced-probe predictions with the 40-probe predictions. 
Results indicate that the number of dynamic probes can be 
reduced to as few as 2 and still retain good accuracy. 
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SYMBOLS 


APD 

a 

2L 

b 


erf (x) 


amplitude probability density 

radius of vortex core (radius at the point 
where maximum velocity occurs) 

mean vortex core size 

radial weighting factor, see table I 

diameter of ring j , see table I 

the error function of x, 

erf (x) = 7 ^/ 0 X e" C d£ 


f 

£' 

*c 

G ^Pt /q 2 


frequency, Hz 

intermediate frequency used to match theoretical 
and analytical PSD functions 

characteristic engine response frequency 

real part of the power spectral density function 
of the total pressure fluctuations 


H 


function describing filter characteristics 


IDC 

IDR 

K 

K A2 » K 9 » X RAD 

K D2 

K inst 

XT 

max 

^max-inst 

x ss 

K 


circumferential distortion factor, (see table I) 

radial distortion factor, (see table I) 

generic distortion factor 

combined, circumferential & radial distortion 
factors respectively, (see table I) 
distortion factor, (see table I) 

generic instantaneous value of K 

generic maximum value of K used in the beta 
distribution function 

most probable maximum value of K^ nst 

generic steady-state value of K 

mean value of Ki nst 
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AK 

k 

m 

N 

n 

OD 

P 

PSD 

p l 

P 2 

P 

Pt 

Ap t 

q 

R Apt /q 2 

R'P 

*v 

r 

rms 

r p 

T 

t 

u 2 

u 


difference between K and K ss , AK = K-K ss 

normalizing fact9r in beta probability density 
function 

exponent in beta probability density function 
total vortex flux, number per second 
exponent in beta probability density function 
outer duct diameter, Isee table II 
probability density function 
power spectral density 
cumulative probability 

probability density function for peak distortion 
factor 

static pressure 
total pressure 

total pressure fluctuation relative to the mean 
dynamic pressure, (q=£pv 2 > 

autocorrelation function of the total pressure 
fluctuations 

compressor tip or duct radius 

+ Zy* , shortest distance from probe to vortex 
radius from vortex centerline 
root mean square 

radius from vortex center line to probe 

period of time over which the maximum instan- 
taneous is to be determined 

time 

mean duct velocity 

x component of velocity induced by the vortex 


iii 


u x 


v 

v 0 

w 

X, Y,Z 

»^p » ^p 

Xy »Xy , Zy 


x,y,z 
x’ ,y* ,z' 


0 

r(n) 


Y 

9 


5 

P 

a 

T 


X component of velocity induced by the vortex 
at the probe 

y component of velocity induced by the vortex 

tangential component of velocity induced by the 
vortex 

z component of velocity induced by the vortex 

fixed coordinate system at the compressor face, 
(figure 9) 

coordinates of a compressor face pressure probe 

location of the x'y’z’ coordinate system rela- 
tive to the X,Y,Z system, (figure 9> 

local coordinate system of the vortex, (figure 97 

coordinate system parallel to the X,Y,Z system 
but with the same origin as the x,y,z system, 
(figure 9) 

vortex orientation angle between the y-axis and 

the x’-y* plane, (figure 9) 

denotes the gamma function of n, 

rCn) = / 0 e^x 11 " 1 dx, r(n+l) s nr(n), 

(figure 3) 

angle between the x’ and x axes, with the x 
axis in the x’-y* plane, (figure 9^ 

circumferential coordinate at the compressor face 

nth central moment (taken about the mean) of a 
statistical distribution 

dummy variable 

fluid density 

standard deviation or rms of a time -dependent 
function 

time delay in the definition of the autocorrela- 
tion function 


IV 


Subscripts 

a vortex core size 

cir circumferential 

1 index of rake of pressure probes 

inst instantaneous value 

j index of ring of pressure probes 

K distortion factor 

max maximum value 

min minimum value 

2 mean flow at the compressor face station 

Ap t total pressure fluctuations 

rad radial 

ss steady-state condition- 

v vortex or tangential velocity- 

<*> evaluated for all frequencies 

Superscripts- 
D mean value 

(~) variable is a function of a number of random 

variables 
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INTRODUCTION 


It has long been appreciated that steady-state inlet 
flow nonuniformities (flow distortion) may have a significant 
effect on the stability of a gas turbine engine and may fur- 
ther induce engine stall (surge). Traditionally, steady- 
state nonuniformities are measured at the inlet/engine inter- 
face (compressor face) by using 4 0 or more steady-state total 
pressure probes. These measurements are combined in such a 
way as to produce a single distortion factor according to 
different definitions (some common ones are listed in 
Table 1) . These distortion factors are then correlated with 
the surge level. 

It was not until recently that attention has been paid 
to the effect of dynamic distortion which reveals the time- 
variant characteristics of the distortion pattern measured 
at the compressor face. This dynamic effect may cause 
engine surge while the steady-state distortion is substan- 
tially below the surge level. Therefore, determination of 
the dynamic effect of inlet flow becomes more important if 
high performance inlets are to be attained. 

Air flow convecting through the inlet duct is always 
unsteady due to flow disturbances generated internal and 
external to the inlet. These disturbances include exter- 
nally generated turbulence. External turbulence is caused 
by atmospheric turbulence, the flow field of the aircraft 
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upstream of the inlet and by external shock waves when in 
supersonic flight. Internally generated disturbulences are 
caused by non-uniformities due to change in the contours 
of the inlet duct and by shock boundary layer interactions. 
Experimentally, this flow unsteadiness is measured at the 
compressor face in terms of total pressure fluctuations by 
the use of 40 (commonly used) or more high response total 
pressure probes (dynamic probes). These unsteady total 
pressure data are analyzed deterministically in terms of 
time-variant distortion factor throughout the data record- 
ing period. The maximum peak distortion factor value of 
the dynamic distortion factor is then screened and used as- 
a design reference. "This method of data analysis is quite 
tedious and expensive. There is a real need to develop 
methods to predict both the peak dynamic distortion factor 
and the corresponding compressor face total pressure pattern. 
These are important for inlet design and data analysis, 
however, they are dependent upon the statistical character- 
istics of the unsteady total pressure data. It is for this 
reason that statistical treatments which utilize the random- 
ness of the total pressure data have been employed to predict 
the peak dynamic distortion factor from limited test data 
(e.g., root mean square, rms , level of total pressure fluc- 
tuations). Improved efficiency and reduction in cost have 
been achieved by this approach. Success of the statistical 
approach is made possible due to the fact that probability 
density distributions of the time-variant distortion factor 
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and the time-variant total pressure data are nearly Gaussian 
(normal distribution) around their mean values (ref. 4) . 

This is true for most of the test cases which do not involve 
inlet flow of strong interaction. Figure 1 shows the physi- 
cal and statistical characteristics of a typical inlet test 
data. A time history of the total pressure fluctuations shown 
in Figure 1(b) is measured by the dynamic probe shown in 
Figure 1(a). Probability distribution of the total pressure 
data is shown to be normal as illustrated in Figure 1(b). 

The time-variant distortion factor resulting from the com- 
bined effect of ail 40 probes is illustrated in Figure 1(c) 
which also shows a normal distribution of the probability 
density of the distortion factor. 

Three major statistical approaches have been developed 
to predict the peak distortion factor and distortion pattern. 
The first one, due to Jacocks, uses Gumbel's extreme value 
theorem to extrapolate the peak distortion value within a 
specified time period from the first short time segment of 
deterministic information (ref. 1). The second approach, 
due to Motycka, et al. , utilizes a random number generator 
for simulating the total pressure fluctuations at each probe. 
The simulated values are based on the measured total pressure 
rms level and the steady-state value. Time-variant distortion 
factors are calculated by using the generated total pressure 
data, from which the peak dynamic distortion factor is de- 
termined (ref. 2). The third approach, developed by Melick, 
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et al., assumes that dynamic effect of the inlet flow is 
totally attributed to the pressure disturbances caused by 
a sequence of convecting vortices with random variables. 

The peak distortion factor is determined from this physical 
model and its statistical properties. The main variables 
of the vortex flow model are identified by matching analyt- 
ical and measured total pressure rms levels (ref. 3)« 

Although all the three methods produce satisfactory 
results, the Melick method is least costly in terms of" 
instrumentation and analysis effort. It can be used online 
while the test is in progress (ref. 5 ). This makes the 
Melick method a very attractive tool for early analysis la 
the inlet design process. 

Also, to minimize costs of dynamic instrumentation and 
overall inlet development, it Is desired to predict the 
dynamic distortion in the early stages of inlet development 
using a minimum number of dynamic probes. Each high 
response total pressure probe is very expensive. This can 
be achieved, in principle, by using the Melick method as 
mentioned in ref. 3 and ref. 5. Melick’ s method has the 
potential to achieve this goal since its peak distortion is 
predicted by the average value of total pressure rms levels. 
Ideally, a single probe would be sufficient to provide the 
same rms level. 

Acceptance of the Melick method has been slow due to 
a lack of understanding by the inlet test and development 
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engineers. Thus, the first objective of the present study 
is to explain and clarify Melick's Technique. The second 
objective. is the development of criteria of dynamic probe 
selection l location) based on the steady-state pressure map. 

A brief review and discussion of the Melick method is 
presented in the first section followed by a series of data 
comparisons using HiMAT inlet test data set (ref, 6). 
Finally, the last section contains the development of the 
criteria for selecting the location of the dynamic probes. 
The results of numerical experiments to validate the probe 
selection criteria are shown. 


THE MELICK METHOD 


Vortex Model 

Physically, total pressure fluctuations of inlet flow 
measured at the compressor face can be interpreted as time- 
variant velocity variations on the compressor face. This 
time-variant velocity variation can again be expressed in 
terms of time- variant vorticity variation. Therefore, it is- 
envisioned that the total pressure fluctuations are caused 
by a sequence of time-variant vortices passing through th<r 
compressor face at a rate of N vortices per second. Each 
vortex would then create total pressure variations at the 
compressor face and generate certain instantaneous (time- 
variant) distortion factors, K inst , at that instant of time 
as illustrated in Figure 2. Figure 2(b) shows that the 
time- variant total pressure, measured by a probe, due to a 
convecting vortex which has a velocity profile as shown ia- 
Figure 2(a>. 

It is assumed that the instantaneous distortion factor, 
Kinst> cau sed by a single vortex is treated as a single dis- 
crete event. Thus, there are NT discrete events during the 
data recording period T. Each vortex has an arbitrary location, 
orientation, and strength. It is further assumed that these 
discrete events are random and independent drawings from a 
probability density distribution described by a Beta distribu- 
tion function (ref. 7). The reason for choosing the Beta dis- 
tribution is because it is the only distribution function 
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that has one bounded positive variable and the ability to 
modify the probability distribution by varying two simple 
parameters. Figure 3 shows the general characteristics of 
the Beta distribution function including the effects of the 
parameters, m and n , on the shape of the distribution. 

The variable of the distribution, x, is shown to be bounded 
within the interval of 0. and 1. The shapes of the dis- 
tribution are skewed to the left in Figure 3 when n^ is 
greater than m . An opposite skew would result if n is 

X X 

less than m^. V/ hen n x and are identical, a Gaussian 
(normal) distribution can be obtained. 

Using the Beta distribution, the probability density 
function of the instantaneous distortion factor, K inst , is given 



where K inst; is rescaled (normalised) by its upper bound, 
K maY , and Kg is a function of m^ and n^ through the gamma 
function as shown in Figure 3. Equation (1) involves three 
unknowns, K^, mg and ng, which require three equations 
for solutions. These three equations can be obtained by 
using statistical moments analysis for P (i.e., the first 
moment, the second moment and the fourth moment, ref. 3). 

The first and second moments of P are related to the mean, 

K, and variance, dg, of the instantaneous distortion factor, 
^inst» respectively. The fourth moment (Kurtosis) of P is 
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assumed to be zero since the distribution of P has been 
shown (ref. 4) to be nearly Gaussian (normal distribution) 
for most experimental data. The parameters, K and d^, can 
be obtained from an analytical expression of K^ nst . and its 
statistical moments. These can be derived from the vortex 
flow model with the main variables, a (mean vortex size) and 
N (vortex flux) being identified by matching analytical and 
measured total pressure rms levels. This will be discussed 
in detail later. Assume, for the moment, that K and d^ 
are already obtained and the distribution of P is known. 
Figure k shows some illustrations of the distributions of 
P with various values of K and d^. Note that the variable 
of P, Kinst, has been normalized by K and d^. That is. 



Equation (la) is the essence of the Melick method, 
from which the peak distortion (i.e. near the right end of 
Figure 4) can be extrapolated by using the statistical 
characteristics of the peak values. Unknowns of this model 
must be related to the inlet flow conditions in order to 
rescale the analytical model close to the experimental 
measurements . 

It is of primary importance that the peak value of 
^inst the NT events be determined from the distribu- 

tion of P. Before this peak value can be found, conditions 
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for the occurrence of the peak value must be determined. The 
peak value implies the largest value of among the 

whole population of NT events. In other words, it is the 
only single event that can exceed a certain distortion 
level. This corresponds essentially to the situation in 
which the deterministic peak distortion factor is found. 
Therefore, the probability of one single event that exceeds 
a certain distortion level needs to be defined. This can be 
described by a cumulative probability of P above a certain 
distortion level. Equation (2) defines this cumulative 
probability of P as P. . Results are illustrated in Figure 
5 for various values of K and d^. Note again that K inst has been 
normalized by K and 6 ^. 



Since it is required that only one event, caused by 
one vortex, can create a peak distortion above a certain 
level, the whole population is therefore required to have 
only one event with probability P^ and (NT-1) events with 
probability (1-P^). This is like the Bernoulli test. The 
overall probability for the condition of only one peak dis- 
tortion above a certain level is therefore described by the 
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Bernoulli (binomial) distribution function given by 
equation (3) (ref. 7) and is illustrated in Figure 6 . 

P 2 - (Y) PlCl-Pl)"- 1 = NT Pi CLP!)— (3) 

From equation (3), the objective peak distortion cor- 
responds to the condition when P 2 is maximum for a fixed NT* 
(i.e., the most probable condition for the peak dynamic dis- 
tortion to occur). This condition is found to be 
when P 2 reaches its maximum. This is obtained by letting 
the first derivative of ? 2 (equation 3 ), with respect to 
P^, equal zero. This can be seen clearly when Figure 4- and 
Figure 6 are combined with the aid of Figure 5 to produce 
Figure 7« For example, when NT=10 and K/o^=4, Figure 6 showa 
that at maximum P 2 , Pj= 0.1 approximately. Using this value 
of P^ in Figure 5, the corresponding value of ( K 1 ^ 5 t-K )/<$£. 
can be obtained and shown in Figure 7. In this way K the 
distribution curves of Figure 6 can be transformed into 
Figure 7 through the aid of Figure 5. In Figure 7, P de- 
scribes the general population of Kj_ nst , and ? 2 describes 
the probability of the peaks of K inst . The peak dynamic 
distortion factor is therefore chosen to be the most probable 
value of Kinst among the population of their peaks, P 2 , for 
given values of NT, R and as shown in Figure 7. Results 
of Figure 7 are generalized and expressed in terms of K, dir 
and NT and are illustrated in Figure 8 . From Figure 8 , the 
most probable maximum value (peak value) of distrction 
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factors can be obtained easily for given values of K, d^ 
and NT. 

In Melick’s original formulation, confidence levels of 
the predicted peak distortion factors were not determined. 
However, with the help of Figure 7, confidence levels can 
be determined easily from the distribution of ?£. Note 
that the variance of ?2 decreases as NT increases. In 
most practical applications, NT and K /d^ are large (~10 S 
and -~10 respectively) such that they fall within the shadded 
band of Fifure 8. 

Matching Model to a Particular Test Case 

Up to this point, prediction of the peak distortion 
factor for a particular case has been accomplished for given 
values of the mean and rms values of the instantaneous 
distortion factor, K and and the number of vortices, NT, 
within data recording period T. It remains to be shown 
how R, d£ and N (vortex flux) can be obtained. From the 
test data, total pressure rms values are measured at the 
compressor face. Therefore, it is intended to derive an 
analytical expression for the total pressure fluctuations 
based on the hypothetical vortex flow model. The main vari- 
ables of the vortex flow model, (mean vortex size (a) and 
vortex flux (N)), are then identified by matching the ana- 
lytical and measured total pressure rms levels. In addition, 
relations between the analytical total pressure fluctuations 
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and the statistical properties of distortion factor, the mean 
value (K) and rms value (d^), are established so that K and 
djt can be determined. 

Development of the vortex flow model originates from 
Melick's hypothesis that the cause of the total pressura 
fluctuations is totally attributed to the convecting vor- 
tices at the rate of N vortices per second. For simplicity, 
a one-dimensional incompressible steady vortex flow solu- 
tion is used for the basis of mathematical development. The 
basic equation for the velocity profile as a function of 
distance from the center of the vortex is given by equation 
(4). (See Figure 2 for definition of terms) 


v e 

°max 




(4) 


where a is designated as vortex size at a radius r=a such 


that 7 q =Vq ms y . 7 g represents the strength of the vortex 
which is the maximum swirling velocity around the axis of 
the vortex. This basic model does not describe the real flow 


situation but it is desirable to keep the selected model as 


simple as possible so that the final expressions will not 
involve too many variables and therefore become overly com- 
plex. Since test data are eventually used to identify the 
main variables, simple expressions do not lose too much in 
accuracy. 

An expression of total pressure fluctuations can be 
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derived from the 1-D vortex model by the use of incompres- 
sible Bernoulli's equation and the spatial relationship 
between the vortex and the dynamic probe as shown in 
Figure 9. The final expression of the total pressure fluc- 
tuation caused by a single arbitrary vortex on an arbitrary 
compressor face probe is given by equation (5). (The devel- 
opment of this equation is shown in the Appendix.) 


^ ' ( U 2. v 8 raax .a.VZ v .VZ p>Y ,S,t) 


u x =v e. 


f Y P- Y V • ^D‘^V 

i- sm y + cos ycos i 

'a a 


Z__Z. 


and r p represents the time-variant distance between the 
probe and the axis of the vortex. 

The subscripts p and v refer to the positions of the 
probe and the vortex respectively. The other variables are 
shown in Figure 9 and are described in the symbol section. 

From this expression [ equation 15H and experimental values of 
total pressure rras levels, three separate branches of develop- 
ment are undertaken to determine the mean vortex size (a), the 
vortex flux N, K and the standard deviation of the distortion 
factor, <?£. 

Determination of Mean Vortex Size 

First, the total pressure rms expression is derived 
from equation (5) through autocorrelation function and power 
spectral density (PSD) function transformations (ref. 9) as 
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shovrn in detail in the Appendix. In order to do these 
developments, probability density functions for the variables 
shown in equation (5) need to be defined to account for the 
effect of fluctuating data. The Beta distribution function is 


used for the main variables, a and V 9 


max* 


and a uniform- 


icons tant) distribution is used for those spatial variables 
describing location and orientations of the vortex (rp,rand 
£) . The final expression of the analytical total pressure 
rms relation is given by equation (6}« 


cr 2 Ap t / , 

q 2 _ (filtered rms level) 2 - C rf (7. 98 f -*■ ) (6) 




a 2 Ap t / q> 


(unfiltered rms level) 


U 2 ‘ 


Equation (6) states that the square of the ratio of filtered 
to unfiltered total pressure rms levels can be expressed by 
an error function of (f / a)/U 2 where is the applied filter 
cut-off frequency and is the mean flow velocity at the 
compressor face. Using this equation, the mean vortex size 
a can be found for given values of f' U2 and the ratio of 
filtered and unfiltered rms levels measured by each dynamic 
probe. It is preferred to choose a certain value of f / so 
that the ratio of the rms level of equation (6) is about 
0.5. This is because equation (6) corresponds to a one- 
point curve-fit. An rms ratio of 0.5 assures the best fit 
over the high and low frequency spectra of the PSD function. 

Equation (6) is a numerical approximation derived from 
equation (5) by assuming that all of the vortices have the 
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same size, which is Che mean value, a, of the distribution. 
Validity of this assumption is examined by comparing the 
exact vortex size distribution given by the Beta function 
and the simplified solutions of the autocorrelation function 
in terms of percent error as function of n a and n a (the 
parameters of Beta distribution function for vortex size a). 
This is shown in Figure 10 along with the range in which the 
experimental data fall. 


Determination of Vortex Flux 
To determine the vortex flux , N, a statistical moments 
analysis (the second and the fourth moments) is applied to the 
amplitude probability density (APD) function of the total 
pressure fluctuations shown in equation (5). Since the APD - 
function, for most experimental data, is shewn to be nearly 
Gaussian (ref. 4), both the first moment (mean) of APD func- 
tion and the 4th moment (Kurtosis) of APD function are 
assumed to be zero. Using the 2nd and 4th moments of APD 
function, one simple expression for N is approximated numer- 


ically by further assuming that nax /^l max ^ ** Equa- 
tion (7) showns the result of this development (see Eqn. 17 


Appx. A) . 


N = 0.2S4 


U 2/*T . C^mH) 


- 0.254 


C^/Rt ) 2 


U2/RT 

(a/Ry) 


(7) 


where denotes the radius of the inlet duct at the com- 
pressor face station. The validity of the simplified 
assumption is also examined and illustrated in Figure 11, 
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along with the range in which typical experimental data 
falls. The implications of the discrepancies between theo- 
retical and experimental results will be discussed later. 

Determination of K and 

Using statistical moments analysis, the mean value and 
the rms value of the instantaneous distortion factor, Y and 
a^, can be expressed analytically by using the expression of 
the total pressure fluctuations (equation (5)). From this, 
X and cr^ can be obtained for given values of the mean vortex 
size a, steady-state distortion factor K sa , and the total 
pressure rms level. 

This is done by substituting equation (5) into the 
appropriate distortion factor equation of Table 1 and then 
substituting these into expressions of the first and second 
statistical moments of the instantaneous distortion factor 
(KinsO as shown in the Appendix. Results of the expres- 
sions of K and are solved numerically in terms of a", K ss 
and total pressure rms level as shown in Figure 12(a) 
through (f) for several distortion factors of Table 1. 

These figures show that K and depend on the value of 1, 
K 5S and the measured total pressure rms level ff Apt/q2- For 
given values of a (found from equation (6)), K ss and rms 
level, Y and can be obtained directly from 12(a) through 
(f) . Using these in equation (7), the value of N is 
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calculated and the peak distortion factor can then be 
obtained from Figure 3. 

Details of the above development are shown in the 
Appendix. 

Contribution of Errors 
in n and m to the Final Results 


Two major simplifying assumptions in the Melick method 
has been examined and shown in Figures 10 Sc 11. One of the 
assumptions is that the vortex size concentrates on its mean 
size, a. The other one is that Vg max / ( V e 3 max )* equals to one. 

There is no way that these assumptions can be justified 
by direct comparison with the test data. However, for engi- 
neering application purposes, it is possible to judge the 
validity of these assumptions approximately by using the 
mean values of the Beta distribution function of the main 
variables, mean vortex size (a) and mean vortex strength 
(Vg max ). Since the mean value of the Beta distribution func- 
tion is a function of its parameters (m and n) only, values 
of the parameters can be resolved from the mean values 
predicted in Melick’ s program by the test data. These 
parameters then can be used in Figure 10 and Figure 11 to 
judge the validity of the two major assumptions. Equations to 
do this are (Sa) and (3b) which give the relationship between 
the mean values of the two main variables and the parameters 



(8a) 

(8b) 
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For values of a/E-^ and ^9 max /U 2 calculated in Melick’ s 
program, relationships of their parameters (m a vs. n a and 
my. vs. n v ) can be obtained from equations (3a) and lob). 

For example, from test data, the calculated values of a/R^ i 3 
ranging from .12 to .18 and the calculated values of V 0 
is ranging from .06 to .10, errors due to the two major assump- 
tions of the Melick method can be approximated as shown by 
the shaded area in Figure 10 and Figure 11. Results show that 
a/R<P is over-estimated by about 15 percent and N is under- 
estimated by about 500 percent. The 500 percent error is due 
to the fact that Vq max/ max ^ assuneb to be 1-0 instead 

of 5.0 as indicated by test data. The error in a/R.p is 
acceptable for engineering applications but the error in N is 
alarming at the first glance. However, the prediction of the 
peak distortion factor is not sensitive to the value of N. 

This can be seen clearly in Figure 8 which has a very small 
slope indicating that large changes in N and T can occur 
without seriously affecting the value of K ma Y _ . ^ . Also, the 

value of dj£ is always an order of magnitude less than K which 
further reduces the effect on K ma Yj1 - nqt . . For most of the test 
data, the value of NT is between 10 4 and 10 s . A 500 percent 
error in N can cause only about 2 percent error in the pre- 
diction of the peak distortion factor. This 2 percent error 
is then negligible. This is why the prediction of the peak 
distortion factor by the Melick method is good even with 
such gross error in N. 
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Effect of Engine Response Characteristics 

Depending on its size and inertia, an engine can only 
respond to distortions below a certain frequency. Since the 
engine compressor does not respond to turbulence above this 
frequency level, pressure measurements should be filtered so 
that high frequency distortions beyond the engine’s response 
capability do not produce meaningless results. That is, high 
frequency components are of no concern because the engine 
can’t respond to them, even if they produce distortions in 
excess of that required to cause a surge. This filtering 
effect can be accomplished by using a third-order Buttbrworth 
filter, which is commonly used to simulate the character- 
istics of present gas turbine engines (ref. 10). The general 
filtering process is shown schematically in Figure 13 which 
shows an input time-dependent function Y(t), as passed through 
a filter H(f), produces an output signal Y’(t) with lower 
magnitude. By properly selecting the filter characteristics, 
this filtering process eliminates the pressure or distortion 
factor fluctuations at frequencies beyond the response charac- 
teristics of the engine. Filtering has the effect of reducing 
the total pressure rms level. The lower the cut-off fre- 

quency the lower the rms level. This is seen in Figure 14 
which shows the ratio of the filtered and unfiltered rms as 
a function of the normalized cut-off frequency, f c a/U^. From 
Figure 14, the ratio of filtered to unfiltered rms level can 
be found for the given value of f c a/U2 to account for the 



20 


engine response characteristics. 


Special Treatment for Some Distortion Factors 


Some distortion factors, such as K A2 and APES, are a 
combination of other indices and must be generated within 
the Melick program. It is assumed that these indices are- 
independent statistically. Prediction of the peak distor- 
tion values of these indices is done by finding the peak 
distortion values of the individual components independently 
and then combined into the indices. This approach is based 
on the assumption that the individual components reach their*- 
peak values at the same time which may not be the case in 
the real situation for two statistically independent indices. 
This assumption assures the prediction of the objective peak 
distortion factors will always be conservative. 

To illustrate this, K^ 2 is defined as: 


k A 2 S 6 + bK RAD 


(9) 


where the b factor is a constant. This means Kg and K^^ 
are independent in equation (9). The peak value of K^2> 
therefore, corresponds to the peak values of Kg and Kj^jj 
simultaneously. The peak distortion values of Kg and Kpjy} 
are determined independently in Melick' s program and then 
combined to obtain the peak value of K^2 using equation 
(9). 
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Melick Mapping Methods 

A very simple method for determining the compressor 
face peak distortion map is included in the Melick program. 
It assumes that the linear vortex of size a is oriented 
along the mean shear line (a straight line) between the 
steady-state high and low pressure regions. The direction 
of rotation of this vortex tends to enhance both the low 
and high pressure regions so that the distortion level can 
be magnified. The peak distortion map is then obtained by 
adjusting the strength of the mean vortex until the predict- 
ed peak value of the distortion factor is reached. 

In reality, the core of the distortion- influence vor- 
tex may not be a straight line. This is why the Melick map- 
ping method is not too accurate in predicting dynamic peak 
distortion patterns. Other methods of adding the vertex 
effects need to be investigated to improve the prediction. 

The Melick method so developed provides a simple and 
low cost way of predicting the peak dynamic distortion. The 
validity of the method is illustrated in the next section. 



DATA COMPARISONS 


In this section, the accuracy of the Melick method in pre- 
dicting peak dynamic distortion factor and pattern will be 
demonstrated by direct comparison with the HiMAT inlet test re- 
sults. In ref. 6, an experimental investigation of a subscale 
HiMAT model with forebody, canard and inlet is described. 

These tests were conducted by NASA in the NASA Lewis S’ x 6' 
supersonic wind tunnel (ref. 6). The HiMAT model has an 
under-fuselage inlet with a high-divergence S-shape duct. At 
the compressor face, there were 40 steady-state total pres- 
sure probes and 40 high response total pressure dynamic 
probes installed on the compressor face station as illus- 
trated in Figure 1. Steady-state total pressure, dynamic 
fluctuating total pressure and filtered and unfiltered rms 
total pressure levels were recorded from the compressor face 
instrumentation during the test. In the present investi- 
gation, data obtained in these tests were analyzed by two 
different methods to determine, by direct comparison, the 
validity of the results predicted by the Melick technique. 

First, the fluctuating total pressure data v/ere digi- 
tized and screened, on distortion factor APRS, through tne 
Dynamic Data Editing and Computing (DYNADEC) System of 
iMarous and Sedlock (ref. 10), from which the measured peak 
dynamic distortion factors and patterns were obtained de- 
terministically. Results obtained from this system are 
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designated as measured results. Second, the filtered and 
unfiltered total pressure rms levels were used in Melick's 
statistical method for predicting the peak dynamic distor- 
tion factors and patterns for the same length of data run 
and are designated as predicted results. Steady-state 
distortion factors and patterns are also included in both 
measured and predicted outputs. 

Because of the high curvature and divergence angle of 
the inlet duct, boundary layer separation is observed. 
Fluctuating total pressure data from the tests shows a down- 
ward spiking characteristics for most of the test cases with 
a clean inlet duct (i.e., without vortex generators). This 
is believed to be caused by intermittent boundary layer 
separation and reattachment (ref. 6). The spiking phenomenon 
of the total pressure data, as illustrated in Figure 15, 
shows random downward spikes bounded approximately by the free 
stream total pressure at the high end and the local static 
pressure at the low end (ref. 6). Because of the spiking 
phenomenon, those cases without vortex generators resulted 
in more than one peak distortion pattern as shown in the 
DYNADEC results. A typical example is shown in Figure 16 (a). 

In order to prevent the inlet performance from being 
jeopardized by potential stall due to the spiking phenomenon, 
a set of vortex generators was installed annularly on the wall 
of the inlet duct down stream of the cowl in order to trip 
the boundary layer and prevent flow separation (ref. 6). 
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Spiking data was not present for most of the test cases with 
vortex generators, indicating good control of the boundary 
layer. 

In the DYNADEC system, peak dynamic distortion of the 
H1MAT data set were screened on the distortion factor APRS 
which is a function of IDC and IDR as given by the follow- 
ing equation. 

APRS - K (IDC) b + X (IDR) (10) 

c r* 

K c is the circumferential sensitivity factor, X r represents 
the radial sensitivity factor and b factor is a function or 
IDR at the tip of the compressor blade (IDR t ^p). This APRS is 
not predicted in Melick' s program. In order to make compari- 
sons possible between measured and predicted peak distortion* 
APRS needed to be computed from the predicted IDC and IDR.. 

It was found the maximum value of APRS occurs at the instant 
when IDC is very close to its maximum and IDR is near its- 
minimum value. This is due to the characteristics of the 
b factor, however, the Melick program does not predict 
minimum instantaneous distortion factor. A study was there- 
fore taken to determine the minimum IDR following Melick r s 
approach. The concept of the prediction of the minimum 
IDR is similar to that of finding the maximum IDR. This is 
illustrated in Figure 17 which shows that the maximum value 
and the minimum value of the predicted instantaneous distor- 
tion factors are symmetric to the mean of the whole popula- 
tion, X. Results were then used for c cnipu ing the peak value 
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of APES. However, note that the Melick program does not 
have the capability to predict peak distortion pattern 
based on PRS. Therefore, the predicted distortion pattern 
are still based on the maximum IDG*. 

The range of the test conditions in the data base are 
listed in Table 2. Data comparisons were made by comparing 
distortion factors and patterns. The first part was the 
comparison of distortion factors. The steady-state values 
obtained from both systems are identical for every case as 
would be expected. Comparisons of the measured and predicted 
peak dynamic distortion factors, APRS, for cases with and 
without vortex generators are illustrated in Figure 18(a) and 
18(b) respectively. Figure 18(a) shows satisfactory agree- 
ment between measured and predicted results. Most of the 
predicted data shown in Figure 18(a) are about 20 percent 
above the measured results, which are on the conservative 
side of the data prediction. This is desirable for an inlet 
design purpose. On the other hand, large discrepancies are 
revealed in those cases without vortex generators, which is 
shown in Figure 18(b). This is because the Melick method 
can not adequately deal with inlet flow when intermittent 
separation induced spiking occurs. An improved method is 
needed to handle the effect of flow separation. An attempt 
was made to use data measured by those probes that were locat- 
ed in the low pressure region (i.e., where high turbulent 
fluctuations in terms of rms level was observed) to make the 
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predictions of the no- vortex-genera tors cases close to the 
measured ones, however, no consistent result was obtained. 

The second part of the data comparison was concentrated 
on the distortion patterns. Again, the steady-state distor- 
tion patterns from both sources were exactly the same. For 
cases without vortex generators, peak distortion patterns 
were not comparable when the spiking phenomenon was observed. 
A typical example is shown in Figure 16. Comparisons of 
the measured and the predicted peak distortion patterns for 
some representative cases with vortex generators (without 
spiking) are shown in Figure 19(a) through 19(n). Although 
the predicted peak distortion patterns are based on maximum 
IDC, instead of APRS, good agreement of most of the compari- 
sons are shown. This may be due to the fact that IDC is the 
dominant term in APRS. 

In some cases, Figure 19(a), 19(d), 19(j), and 19 (m>- 
for instance, discrepancies between measured and predicted 
patterns are the most apparent. It may be possible to im- 
prove the pattern prediction by changing the axis of the 
main influence vortex from a straight line to a more realis- 
tic coutour line as discussed previously. More effort is 
required to improve the prediction. 


METHOD OF PROBE SELECTION 


Development of Probe Selection Criteria 

In principle the Melick method requires only one dynamic 
total pressure rms measurement providing it is located at a 
point that represents an average rms level for the whole com- 
pressor face ti.e., 40-probe average). If the number of 
probes required to provide this average rms level can be mini- 
mized, the cost of test instrumentation and perhaps the flow 
blockage of the inlet can be reduced greatly. 

If it is assumed that dynamic effects of the inlet flow- 
are somehow related to the steady-state total pressure pattern, 
a basic rationale of dynamic probe selection can be established 
based on the steady-state pattern. This basic assumption 
provides a foundation for this investigation. 

Since the peak dynamic distortion factor from the Melick 
method is based on the average values of total pressure rms 
level and mean vortex size, a, at the compressor face, an 
examination of the total pressure rms levels at the compressor 
face should reveal any similarity between the steady-state 
and dynamic characteristics of the inlet. The relationship 
between the steady-state total pressure and the rms levels 
are examined by plotting radial variations of both values as 
shown in Figure 20(a) through 20(f). It has been found that 
the rms levels have an inverse relationship to the steady- 
state values. That is, when the steady-state total pressures 
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are high, the dynamic content of the flow is low and vice 
versa. It is further noted that when the rms value crosses 
the average value (rms/rms^ = 1) the steady-state total 
pressure also crosses the mean value (P t /P c+<J = -) • This- 
observation suggests the possibility of establishing a cri- 
terion for selecting apriori the location of dynamic probes 
that will be close to the mean rms level of the duct. Never- 
theless, the rms levels do not represent the entire dynamic 
effect because another variable, mean vortex size a, also- 
has an important effect on dynamic behavior. 

Two approaches were taken to establish probe selection 
criteria. The first approach was done by searching for com- 
binations of dynamic-probe subsets which produce average values 
of rms level and mean vortex size equivalent to the 40-probe 
average values. This approach was found to be too complex 
because of the difficulty in handling two variables in conjunc- 
tion with a multiplicity of combinations of the dynamic probes. 
A second approach was employed by examining local con- 
tributions of dynamic effect to the peak dynamic distortion 
factor directly. In order to do this, the Melick program 
was modified so that the peak distortion factor could be 
predicted by using dynamic data measured by each individual 
probe. Results were tnen compared with the peak distortion 
factor predicted by all 4-0 probes ^uO-probe value). Radial 
variations of these local contributions to peak distortion 
factor IDC are also plotted on Figure 20. Similarities 
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between the trend of the ras levels and the distortion fac- 
tor variations are remarkable as seen in Figure 20. The 
single-probe distortion factor variations can be interpreted 
as variations of overall dynamic effect (including the effect 
of rms level and mean vortex size a) over the compressor face. 

Only IDC was used throughout the development since IDC 
is the dominant term in computing APRS which is the basis 
for comparison with the DYNADEC results. Results of the 
present investigation are also applicable to other distortion 
factors since the influence of the dynamic component to every 
distortion factor is the same in the Melick method. 

In order to find the key probes, the line of perfect 
agreement between the single-probe value (of peak distortion 
factor) and the 40-probe value were plotted on the steady- 
state pressure map which is the basis of dynamic probe selec- 
tion. Figure 21(a) through 21(1) show this comparison for a 
number of operating conditions. Note that only the steady- 
state average pressure contour and the dynamic average line 
were plotted in Figure 21. This was done for most of the 
HiMAT test cases with and without vortex generators and 
covered the whole range of the test Mach number, angles of 
attack and sideslip. The results show that the lines of 
perfect agreement (dynamic average lines) coincide fairly 
well with the steady-state average total pressure lines 
except in regions of high and low pressure. In other words, 
the average turbulence line, i.e., the line of perfect 
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agreement, is in general in the vicinity of the average 
steady-state total pressure region. This implies that 
dynamic probes selected close to the steady-state average 
line would give the prediction of a peak distortion factor 
very close to that of the 40-probe prediction. It would be 
preferred that the predicted value always be on the conser- 
vative side, i.e., slightly higher than the 40-probe value. 
This occurs for dynamic probes located on or outboard 
(toward the low pressure region) of the steady-state average 
line. On the other hand, agreement between steady-state and 
dynamic average lines do not hold consistently for regions 
of particular high and low pressure on the steady-state 
pressure map. This can be examined in Figure 19, Figure 20 
and Figure 21. Using the results of the above analysis, 
general criteria for the selection of the dynamic probe 
locations can be drawn. 

General Criteria: 

(1) The dynamic total pressure probes that are on or 
just outside (larger radius) of the steady-state 
average total pressure line should be selected. 

(2) Probes selected should avoid very high and low 
steady-state pressure regions. 

Accuracy of the Criteria 

The accuracy of the above criteria has been tested by 
using the HiMAT test data and satisfactory results were 
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obtained for engineering application purposes. One example 
of dynamic probe selection is shown in Figure 22 for test 
point 194S ( M=0 . 9 , <x=7° and £=0°). The steady-state average 
pressure contour is shown in Figure 20(c) and the associated 
selected dynamic probes are also illustrated according to the 
general criteria of dynamic probe selection. The selected 
circumferential S probes are: nos. 1, 7, 12, 20, 24, 30, 32 
and 37. The peak distortion factor (IDC) predicted by these & 
probes is 0.0470 which is +4.73% higher than that predicted 
by 40 probes. When only four probes out of the selected S 
probes are used (i.e., nos. 1, 12, 24, and 32 ) the predicted 
peak IDC is +1.12% higher than the 40-prooe value. If the 
number is further reduced to 2 probes (nos. 12 and 32 ), the 
predicted peak IDC is +5.10% higher than the 40-probe value. 
Following the same procedure, results of some typical cases 
with vortex generators are tabulated in Table 3. Results 
show that the accuracy for using eight probes circumferen- 
tially located near the mean steady-state pressure is within 
5 percent, whereas selecting the wrong probe could cause 
errors of ±20 percent as illustrated by the single-probe 
worst cases. Table 3 also presents satisfactory results for 
using 2 and 4 probes. The results in Table 3 also fall in 
the histograms of Figure 23. Figure 23 shows a histogram of 
all of the probe selection data analyzed with and without 
vortex generators. With vortex generators, which are repre- 
sentative of conventional low divergence ducts with moderate 
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bends, the results are generally conservative compared to 
the 40-probe result. This is especially true of the 2-probe 
results. However, without vortex generators (a typical high 
divergence duct with bends or obstructions) the results are 
still conservative but the 2-pro’oe result is the least con- 
servative . 

Pragmatically, the locations of the dynamic probes are 
preferred to be fixed during the test so that a sequence of 
test runs can be conducted without stopping to relocate 
instrumentation on the model. Steady-state total pressure 
contours can be obtained from prior developmental test data, 
an analytical prediction or the intuition of an experienced 
inlet/propulsion engineer. In some cases it may be possible 
to run some initial tests and then install the dynamic 
instrumentation later in the wind tunnel test sequence. 

Once the general shapes of the steady-state pressure maps 
are obtained, fixed locations of the dynamic probes can then 
be determined by applying the general criteria. 

The above method of utilizing fixed dynamic probes has 
been examined by using the HiMAT data set ( with-vortex- 
generators cases). Two fixed probe locations were determined 
from the general shapes of the steady-state average pressure 
contours shown in Figure 21. Probes no. 12 and 32 are the 
most promising candidates because they avoid the extreme 
high and low pressure regions and they usually lie on the 
low pressure side of the generalized steady-state average 
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pressure contour. Note that these lie on the opposite side of 
the approximate line of symmetry connecting the regions of 
high and low pressures. Generally, this cross-line-of- 
s ymme try combination performs surprisingly well since these 
two probes usually satisfy the general criteria (1) and (2), 
especially that of avoiding high and low pressure regions. 
Results of the fixed 2-probe cases are also included in 
Table 3 and Figure 23 and are comparable with the other 
results. 

In order to make sure that the general criteria are 
also valid for other inlet test data, numerical experiments, 
of the general criteria were conducted by using 3 data sets 
of Melick’s test cases for an internal compression super- 
sonic inlet. Results and test conditions of this study are 
listed in Table 4. Although Melick’s test cases do not pro- 
vide all the 40-probe data (only 14 dynamic probes were used ) , 
results of the selected 3 or 4 probes show satisfactory 
agreement indicating that the criteria are also good for 
other types of inlets. However, other sources of data 
should be examined before this criteria can be fully vali- 
dated . 

The general criteria provides a way for selecting cor- 
rect locations of the dynamic total pressure probes by simply 
looking at the steady-state total pressure map. Probe reduc- 
tion becomes possible without significant loss in accuracy. 
The number of probes can be reduced to as few as 2 instead of 
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40. Thus, a great deal of savings in instrumentation and 
analysis costs can be achieved. 


CONCLUSIONS 


A simplified description of the Melick method in pre- 
dicting peak dynamic distortion at the compressor face of a 
gas turbine engine has been accomplished. The main variables 
of the assumed inlet flow vortex model have been identified 
by matching analytical and measured total pressure rms 
levels. Statistical properties of the instantaneous distor- 
tion factor are then related to the main variables, from 
which the most probable peak value of the instantaneous dis- 
tortion factor can be extrapolated. The accuracy of the 
Melick method has been demonstrated by using a HiMAT inlet 
test data set and comparing it to the measured results of 
DYNADEC . 

The criteria for dynamic probe selection have been 
established through a sequence of comparisons between a 
steady-state total pressure map and turbulence behavior at 
the compressor face. Instead of using 40 dynamic probes, 2 
probes will suffice in the early stages of inlet testing 
with errors of approximately 5 percent compared to the 
40-probe results. 

The Melick method cannot properly deal with turbulence 
caused by severe intermittent flow separation and does not 
predict better peak distortion patterns than Motycka’s method 
does. However, the pattern prediction can be improved with 
suitable modifications in the computer program. 
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Another analytical flow model and/or modified statis- 
tics are required to handle the effect of intermittent flow 
separation. On the other hand, other methods of construct- 
ing the peak distortion, pattern from the effect of mean 
vortex flow field may be proposed for better prediction of 
the peak dynamic distortion map. 

A logical extension of the current work would be to 
predict the rms dynamic total pressure fluctuations directly 
from the steady-state pressure maps so that dynamic total pres- 
sure measurements would not be needed. This would allow the 
Melick technique to be used as a design tool to predict the 
dynamic distortion before the first developmental tests 


were run. 


REFERENCES 


1. Jacocks, J. L., et al.: STATISTICAL PREDICTION OF MAXI- 

MUM TIME VARIANT INLET DISTORTION LEVELS. Arnold Engi- 
neering Development Center, (NTIS) AD/A- 004- 104 , 

January 1975. 

2. Motycka, D. L. : DETERMINATION OF MAXIMUM EXPECTED 

INSTANTANEOUS DISTORTION PATTERNS FROM STATISTICAL PROP- 
ERTIES OF INLET PRESSURE DATA. AIAA Paper 76-705, 

July 1976, 6 pp-. 

3. Melick, H. C. , Jr., Ybarra, A. H. and Beneze, D. P.: 
ESTIMATING MAXIMUM INSTANTANEOUS DISTORTION FROM INLET 
TOTAL PRESSURE RMS MEASUREMENTS. AIAA Paper 78-970, 

July 1978, 15 pp . 

4. Stevens, C. H. , Oliphant, R. C. and Spong , E. D.: EVALU- 

ATION OF A STATISTICAL METHOD FOR DETERMINING PEAK INLET 
FLOW DISTORTION USING F-1S AND F-18 DATA. AIAA Paper 
80-1109, July 1980, 11 pp. 

5. Sanders, M. E. : AN EVALUATION OF STATISTICAL METHODS 

FOR THE PREDICTION OF MAXIMUM TIME- VARIANT INLET TOTAL 
PRESSURE DISTORTION. AIAA Paper 80-1110, July 1980, 

11 pp. 

6. Neumann, H. E., Povinelli, L. A. and Coltrin, R. E.: AN 

ANALYTICAL AND EXPERIMENTAL STUDY OF A SHORT S- SHAPED 
SUBSONIC DIFFUSER OF A SUPERSONIC INLET. AIAA Paper 
80-0386, January 1980, 12 pp . 

7. Milks, Samuel S.: Mathematical Statistics. John Wiley 

6 Sons, Inc . , 196 3. 

8. Marous , J. J. and Sedlock, D.: DYNAMIC DATA EDITING AND 

COMPUTING SYSTEM (DYNADEC) . Proceedings of the Air Force 
System Command Science and Engineering Symposium, AFSC- 
TR-003. Vol . 1, October 1973. 

9. Bendat, J.S. and Piersol, A. G.: MEASUREMENT AND ANALY- 

SIS OF RANDOM DATA. John Wiley 6 Sons, Inc., 1968. 

10. Marous, J. J. and Sedlock, D. : AN ANALOG EDITING SYSTEM 

FOR INLET DYNAMIC FLOW DISTORTION, DYNADEC- -PAST, PRES- 
ENT AND FUTURE. AIAA Paper 80-1108, July 1980, 10 pp . 

11. Discus, J. H.: FORTRAN PROGRAM TO GENERATE ENGINE INLET 

FLOW CONTOUR MAPS AND DISTORTION PARAMETER. NASA TM X- 
2967, 1974. 


37 


33 


12. Melick , H. C. , Jr.: ANALYSIS OF INLET FLOW DISTORTION 

AND TURBULENCE EFFECTS ON COMPRESSOR STABILITY. NASA 
CR- 114577 , 1973. 



APPENDIX 


Vortex Flow Model 

As described in the main text the total pressure fluctua- 
tions of inlet flow is attributed to the existence of random 
vortices convecting down the inlet duct. For sim- 

plicity, the total pressure fluctuations caused by a one- 
dimensional vortex flow is considered. Then the effect of 
many random one-dimensional vortices will be considered by 
integrating over all possible values of the variables. 

The tangential velocity of one-dimensional, steady, 
incompressible vortex is given as: 

V Q =* Vq (r/a) e -1 / 2 [( r / a ) 2 - l-] 

° a ZL3LX ( 1 ) 

where Vq is the tangential velocity at radius r and V g max 

represents the maximum Vq at radius a. The size and the 

strength of the vortex are defined to be "a" and V ft 

w max 

In reality, the effect of vortex decay due to viscous 
dissipation results in decreased strength and increased size. 
However, in the present analysis the vortex size and strength 
are assumed to be constant because the time required for the 
vortex to travel through the inlet duct is very short and the 
vortex decay rate is rather small. 

The total pressure fluctuations caused by the vortex 
can be superimposed on the steady-state total pressure com- 
ponent to form the time-variant actual total pressure as 
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indicated by the following equation. 


P t + AP t 


where, by the incompressible Bernoulli’s equation, 


( 2 ) 


P T - P s + *P(n 2 + u x ) 
P t “ P s + if 0 ! 


(3a) 

(3b) 


where p^, represents the time-variant actual total pressure, 
P t is the steady-state total pressure, p t is the total 
pressure fluctuations, p s . is the local static pressure, U 2 
is the mean stream velocity at the compressor face station 
and U x is the perturbed velocity caused by the vortex. 

By substituting (3a) and (3b) into (2) and dividing 
through by the mean dynamic pressure q 2 =£pU 2 . An expression 
for the total pressure fluctuations is obtained: 


APt. AP t 20^ 0* 2 

iPuf 1 2 = “ 2 V 


(4) 


Since U x is assumed to be much less than U 2 , the second term 
on the right hand side of (4) can be neglected. 


AP t^ 2U x 


<12 


(5) 


Before equation (5) can be used for further development, an 
expression of U x needs to be derived. This can be done 
easily by relating U x to Vq of the one-dimensional vortex 
through the geometric relationship between the vortex and 
the pressure probe location. This is illustrated in 
Figure 9. 
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For simplicity, all vortices are considered to be effec- 
tive only when they are crossing the measuring plane and the 
pressure probe is located at the center of the duct. After 
a prolonged mathematical development involving geometric 
transformations, the following expression for U x is obtained. 


U 


81 


max 


r-v* 


Zainr+?E^oareosp)e~* (6a) 


where 


l^) 2 = 


-J^cosYcos 0-^2 sin 


) + (_zI_£.)_ (l_|LcosYcosp 

f 


(6b) 


(Yp,Zp) represents the probe location. 

(Y V ,Z V ) is a point on votex axis which has the shortest 
distance to the probe 
-tt/2 < y <. tt/2 


- 7T <_ 8 <. T 

With equations (6a) and (6b) , equation (5) can be 
written in function form: 



flU-jV© > a > »^n > ^ 
w max y * 


( 7 ) 


Next, the vortex model will be used to develop an 
expression for the power spectral density (PSD) function so 
that test data can be compared. To do this an autocorrelation 
function will be formulated by the use of equation (7). 
Finally, randomness of the variables will be taken into 
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account in the formulation by employing the probability 
density functions for the variables from the physical consid- 
eration of the vortex size, the vortex strength* the vortex 
orientation and the vortex location* etc. 

The autocorrelation function of the total pressure fluc- 
tuations due to a single vortex is defined as: 


V/q (T,a ’' r 9 aax' R Y> T > S) =i< T > 

t /H 2 jTi 


t+Tjdt 

q 2 ^2 




^ A 

where R and Ap t represent the functions that involve the 
random variables (a, v e max » p v > Y » P) and ^v = ^ Y v“^p^ ? + 
(Z v -Z p ) 2 . 

After the probability density function for the variables 
are used in equation (8) and assuming that there are N 
vortices per second of vortex flux traveling down the inlet 
duct, a complete autocorrelation function of the total 
pressure fluctuations can be obtained as shown in equation (9K 


^Pt/q 2 (T > N?U)P(Ve max ) 

all possible values 

P(R V )P(T)P(P) dad7 0max dR v drdfi (9) 

where P(a), p ( v 9 max ) » p (Ry) . p (T) and P(p) represent the 
probability density functions for the vortex size a, the 
vortex strength the vortex location R v and the vortex 


orientation Y Sc. £ respectively. Here, it is assumed that all 
the variables are independent. 

The probability density functions for fly, y and £ are 
assumed to be distributed uniformly (i.e. constant). Beta 
distribution function is employed to describe both a and 
V 0 max . The beta distribution function is given as: 

PU) = K ll)°u-£) n , for 0 <5 <1 

5 and m,n £ 0 

With these assumptions, equation (9) can be normalized 
and solved numerically in one of two ways. The first is a 
complete numerical integration which shows that the results 
depend on the parameters m a and n a . The second is by as- 
suming a delta function for P(a) instead of using the beta 
distribution function. The latter results in a single curve 
which is independent of m a and n a , and is dependent only on the 
mean size of the vortices, 5. For simplicity, the results of 
the second integration will be used for further development 
in the analysis. 

With this numerical solution of the autocorrelation 
function, a PSD function can be obtained easily through a 
Fourier transformation. Since only the real part of the 
transformation is comparable to test data, the real part of 
the PSD function is shown in the following equation. 




( 11 ) 
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Equation (11) can be normalized and solved numerically. 
The numerical results are approximated by an exponential 
function and the final expression is given by: 

°APt/q? „ , - 63 . 6 US/U ,} 2 

-5 - 9 . 0 (a/U 2 )e 41 ( 12 ) 

a APt/q a 

Equation (12) is then integrated over part of the fre- 
quency spectrum. The left hand side of equation (12) becomes 
a ratio of the filtered and unfiltered mean square value of 
the total pressure fluctuations. The integration of the right 
hand side of equation (12) results in an error function which 
is described by equations (13a) and (13b>. 


/ G A Pt /q 
\ -= ~d£ 

0 A V<l2 



-63 , 6 ( fa/U 2 )~ 


df 


(13a) 


or, 


erf(7.9#, ll3b 

°ip C /q 2 

Equation (13b) is the key equation for finding the main 
variable a. The ratio on the left is measured in the inlet 
test for a selected filter Cut-off frequency f'. The mean 
flow velocity at the compressor face, U 2 , is also measured 
in the test. 5 can be solved by using equation (13b) with 
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the measured data at each probe. This is done on a probe by 
probe basis and then an average value is taken. 

APD Function Analysis for 
Finding the Vortex Flux M 

Having found a, the vortex flux N remains to be identi- 
fied. To find solution for N, statistical moments of the 
amplitude probability density (APD) function are employed. 
Since experience with test data indicates that the ampli- 
tude distribution of the total pressure fluctuations is 
nearly Gaussian, its kurtosis (the 4th moment) is assumed to 
be zero. It will be shown that the second and the fourth 
moment of the APD function can be used for solving N and 

^Qmax* 

First, the nth moment of the APD function is defined 

as: 

“a “ * J J J J j J ^f] n PCa ) P(V 9 max ,P, V 

all possible values 

P(T) P(P) dtdadVp dR T ,dTd0 mm 

®max v v 

Relations between p n and the statistical properties of inlet 
fluctuating pressure data can be described by equations (15). 
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mean * p]_ * 0 

(15a) 

mean square =■ =* 0 2 

(15b) 

skewness = ^ 3 /(^ 2 ^*^ 

( 15c) 

kurtosis - p.^/u | - 3*0 

( 15 d) 


Equation (14) is solved numerically and normalized for 
the second and the fourth moment of the APD function, & 
p^. Results are approximated as follows. 


£2 

(nr t /u 2 )(v 0jmx /u 2 )* 


^V u 2 Hv 6max / u 2 ) 4 


i^o/KS/r ^ 2 

2 * 75 if u/R T )2 


( 16a) 


(16b) 


Using equations (16a), (16b) and assuming that the kurtosis 
is zero, the following expression is obtained. 


N = 0.254 


U 2 /R.p 

U/Rji ) 2 




■max 


( 7 Z } 2 

' v 0 max ; 


(17) 


The last term of the above equation is assumed to be 
1.0, which is accurate to within 20 percent when m v > 4 + .4n y . 
Therefore, equation (17) becomes: 


U 2 /Rt 

(a/R-j ) 2 


N = 0.254 


(18) 



Equation (18) is the solution for vortex flux N since 
3 is already solved in equation (13b) and U2 and H-p are known 
values. ^Sfnax is then obtained by an algebraic solution of 
equation l 16a). 


Distortion Factor Determination 

Up to this point, the main variables of this vortex 
flow model have been identified using the key equations and 
the experimental data. The next part of the development is 
to establish relationship between these variables and the 
dynamic distortion factor of the inlet. The main variable* 
of the vortex model (I and N) are useful for predicting the 
statistical properties (mean and standard deviation) of the 
distortion factor. From this the extreme (peak) value of the 
dynamic distortion can be determined by extrapolation. 


Relation Between Total Pressure 
Fluctuations and the Distortion Factors 

Typical distortion factors for the inlet flow are 
defined in table I. These may be used to discribed the var- 
iation of pressures at any instant of time. The time-variant 
total pressures can be characterized in terms of the vortex 
flow model variables, a and N, decribed previously. 
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Following the same analysis as for the APD function, the mean 
and the standard deviation of the instantaneous distortion 
factors, K and d^, are found to be functions of a, the steady- 
state distortion factor (K ] and the total pressure rms 
level. These parameters are already known at this point so 
that K and d^ can be obtained explicitly. See Figures 12(a) 
through (f) for illustrations. 

After K and d^ are obtained, probability distribution 
function for the population of the instantaneous distortion 
factors is defined from which the most-probable peak distor- 
tion factor can be picked. The peak value is found to be a 
function of K, d^ , N and the data record length T, which are 
again known values from the previous analysis. 

Now, relations between the total pressure contour 
measured at the compressor face and the distortion factors 
are linked by the definitions listed in table I. Since the 
total pressure and the distortion factors are time-variant, 
it is more appropriate to used the instantaneous notation. 
Instantaneous distortion factor with the random variables 
involved is denoted as K^ nst which is a function of the in- 
stantaneous pressure contour by the definitions of table I 
as indicated in equation (19). 

®inst - - f <Ps3 + ^t> " 1191 
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AP t in equation (19) is obtained from equation (7) and 
the steady-state total pressure data, P ss > are assumed for 
several typical distortion patterns. Statistical moment 

A 

analysis for K^ nst . is taken so as to relate the mean value 
K and the standard deviation to the main variable a. The 

A 

nth moment of K^ nst is defined in equation (20). 

* ■ -J» inst^ ^ j^pumve^mRvmrme) 

all possible values 

dtdadV 0ma df^dTdS (201 

The probability density functions used in the analysis remain 
the same as those used in the APD function moments analysis. 
Only the first two moments of equation (20) are needed for 
finding K and d^ through the aid of the following relations. 

AK-I_-K SS (21a) 

<T k - U 2 - K 2 ) i (21b) 

The results of this analysis are solved and approximated 
numerically, and are shown in Figures 12(a) through (f) for 
different distortion factors listed in table I. It is shown 
that d £ and AK are functions of a and K ss . For some of the 
distortion factors, the results are functions of a only. After 
a and K ss are fed into these figures, and AK can be de- 
termined. Then K is calculated from equation (21a) for the 
given value of K 3g . 
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Table 1. Definition of distortion factors 
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Table 4. Numerical experiments of the general criteria of 
dynamic probe selection for three Helick's test 
cases of an supersonic internal compression inlet 


Case no. 

Test condition 

Totally 14 dynamic 
probes 

Available selected 
probes 

M 

a 

8 

peak K ^2 

worst 

single 

probe 

(%) 

probe no. 

Error 
based on ; 
14-probe 
value (%)r ■ 

1 

2.5 

5 

■ 

1.1638 

+35.81 

-15.65 

18,21,22, 

38 

-Z.5& 

2 

1.6 

5 

1 

1.1636 

+16.58 

-19.13 

4, 17,21, 
38 

-2.97 

3 

2.5 

5 

H 

0.7438 

+25.67 

-14.91 

4, 20,24 

+4.94 
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(a) HiMat Test Model: 0* 



(c) 


Distortion Factor(resulting from a combination of all probes) : 



Figure 1. Illustration of a typical inlet test model and the physical and 
statistical total pressure data and distortion factor 
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P(x) 



Figure 3. Illustration of the characteristics of Beta distribution 
function 



Figure 4. Probability density function of the instantaneous distortion 
factor, K-; n q r . Note that the results have been normalized 
by K and 









Figure 8. Most probable maximum(peak) instantaneous distortion factor 
as a function of NT, K and 
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Figure 12. Continued 
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Time log(f) Time 


Input Signal Filter Output Signal 

Figure 13. Illustration of the filtering effect for engine respond 
characteristics 



Figure 14. Effect of filter cut-off frequency on distortion factor ms 
level 





Time-sec 


Figure 15. schemetical illustration of total pressure spiking 
phenomenon 


(a) Measured distortion pattem(DYNADECl 



measured peak no. 1 measured peak no. 2 


(b) Predicted distortion pattem(Melick' s) 

Data P t.: 1574 
Mach No . : 0.9 
Angle of Attack: 15 deg 
Sideslip Angle: 0 deg 

Figure 16. Typical measured and pridicced distortion pattern for 
data with soiking phenomenon 
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increase NT 


P or P„ 


' P' ,NT=10Q 


•p,K/a K =4 


■p 2> nt=ioo 
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increase NT 


-most-prob-min 
value for NT=100 


K -K 
inst 


-most-prob-max value 
for NT=100 


Note: V' represents the probability density function 
for the minimum peak condition 

Figure 17. Illustration of the method of finding the minimum peak 
distortion value 


line of perfect 
agreement- 


line of perfect 
agreement • 


APRS 

(predicted) 


.03 .C 

APRS (measured) 

(a) With Vortex Generators 


.03 .06 

APRS (measured) 

(b) Without Vortex Generators 


Figure 18. Comparisons of predicted and measured peak distortion 
value , APRS 
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Data Point: 1965' 

Mach no. = 0.9 

Angle of Attack = -10 deg 

Sideslip Angle = -5 deg 
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Figure 19. Continued 



Steady-state 


Measured (DYNADEC) 
Peak 



Data Point: 1948 
Mach no. = 0.9 
Angle of Attack = 7 deg 
Sideslip Angle = 0 deg 
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Figure 19. Continued 
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1-probe worst case 


1-probe worst case 


(a) With Vortex Generators 


(b) W/0 Vortex Generators 


Figure 23. Accuracy of dynamic probes reduction (using the general criteria 
of dynamic probe selection) based on the 40-probe predicted 
peak distortion values 


PART II 

INLET FLOW DYNAMIC DISTORTION 
PREDICTION — WITHOUT RMS MEASUREMENTS 



SUMMARY 


A two-equation turbulence model, k-e model, is used in 
the prediction of inlet flow dynamic distortion of jet 
aircraft based on steady state total pressure measurements 
only. This turbulence model is solved at the compressor 
face station by using a finite difference scheme. Total 
pressure rms level of the inlet flow is predicted by the 
turbulence model. The Melick statistical method is then 
employed to estimate the peak dynamic distortion based on 
the analytically predicted total pressure rms level. 

A statistical method is developed for the estimation of 
boundary conditions for the turbulent equations. 

In order to solve the set of nonlinear, coupled 
turbulent equations, an implicit formulation is utilized 
such that a set of discretized finite difference equations 
can be arranged into tridiagonal matrix equations for 
efficient numerical iterations. 

Results of the prediction compare well with 
experimental measurements of subsonic, transonic and 
supersonic inlets under various flight conditions. The 
present method can be used in the preliminary inlet design 
phases to reduce the design costs. 
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1. LIST OF SYMBOLS 


Symbols 




G 

IDC 



max 


Iog 10 

M 

o 

ms 


N 

N 


S 

r 


Discription 
Tridiagonal matrix 

Specific heat at constant pressure 

Correlation factor between boundary rnvs: 
level and boundary total pressure loss 

Empirical constant of turbulence model 

Empirical constant of turbulence model 

Empirical constant of turbulence model 

Empirical constant for Kolmogorov-Prand tl 
expression 

Empirical constant for the expression of e 

Relative percent error of each numerical 
iteration 

Column matrix 

Relaxation factor of the numerical scheme 
Column matrix 

Circumferential distortion factor, 

(see Table l) 

Combined distortion factor, (see Table 1) 

Turbulent kinetic energy (ftVsec 1 ) 

Turbulence mixing length (ft) 

Logarithmic function with base 10 

Free stream Mach number 

Mean square value of time variant total 
pressure (psi 1 ) 

Grid size parameter, (number of divisions 
between locations of total pressure probes) 

Number of rings of total pressure probes 
iii 


N 

P 

Number of rakes of total pressure 

probes 

P 

Pressure, in pound per square inch 

(psi) 

P 

Nondimensional pressure 


R 

Gas constant, =1545.33 f t- lbf /Ibmole- 3 R 

r 

Radial coordinate (ft) or recovery 
factor of total pressure probe 


rms 

Root mean square of time variant total 
pressure (psi) 

rms 

Compressor-face-average rms level, 



= /rTN ’ £ j) 2 

W r p l j ,J 


T 

Air temperature (°R) 


TN 

Total number of non-boundary grid 

points- 

CM 

3 

Compressor face average mean flow 
(ft/sec) 

velocity 

u 

Axial component of flow velocity (ft/sec) 

ip* 

Mean square value of axial velocity 
fluctuations (ft 2 /sec z ) 

V 

Radial component of flow velocity 

(ft/sec) 

w 

Circumferential component of flow 
(ft/sec) 

velocity 

X 

Column matrix 


X 

Axial coordinate (ft) 



Greek 

a 

0 

•y 

e 


Angle of attack (degree) 

Sideslip angle (degree) 

Specific heat ratio 

Turbulent kinetic energy dissipation rate 
(ftVsec 1 ) 


IV 



e 


p 


v 



Circumferential coordinate (radian) 
Absolute viscosity (slugs/f t-sec) 
Kinematic viscosity, = p/p (ft 2 /sec) 
Turbulent viscosity, = C^k"/e (ft 2 /sec) 
Air density (slugs/ft 3 ) 

Standard deviation 

Empirical constant of turbulence model 
Empirical constant of turbulence model 
Arbitrary variable 


subscripts 

0 

1 

2 

b 

i 

j 

m 

P 

s 

ss 

t 


Free stream condition 
Inlet throat flow condition 
Compressor face flow condition 
Boundary value, (i.e. near the wall) 
Radial index or tensor index 
Circumferential index or tensor index 
Experimentally measured result 
Analytically predicted result 
Static condition 
Steady state value 
Total or stagnant condition 


v 


Superscripts 

(n) 


(over bar) 


Indicate condition after the nth numerical 
iteration 

Indicate solution of equations (22) & (23) 

Time variant component 

Time averaged value or compressor face 
average value 


vi 


2. INTRODUCTION 


Distortion of inlet flow plays an important role in 
determining the stable operating region of turbojet or 
turbofan engines, especially for highly maneuverable 
military aircraft. Time variant or dynamic distortion 
heightens this problem by creating peaks superimposed on top 
of the steady state distortion as illustrated in the 
adjacent figure. Thus, the occurrence of stall surge or 
other engine instabilities can not be readily predicted. 

Traditionally, time variant total pressures of the 
inlet flow measured in the test are analyzed by several 
methods. These methods can be classified into two 
categories . 

In the first category, time variant total pressures 
measured by high response total pressure probes over the 
compressor face are recorded in the test for a specified 
period of time. Then, these data are screened in the 
Dynamic Data Editing and Computing System (DYNADEC) 
deterministically based on a certain distortion factor 
(ref.l), (definitions of some distortion factors are given 
in Table 1). From the screening procedure, a maximum 
instantaneous peak distortion pattern corresponding to a 
maximum instantaneous peak distortion factor is determined 
(refs. 1,2). This is the most accurate method, but it is 
also the most costly and time consuming method. It will be 
the basis for data comparisons in this work. 
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ILLUSTRATION OF SOME FEATURES OF THE TIME VARIANT TOTAL 


PRESSURES AND DYNAMIC DISTORTION: 


Compressor Face 
Instrumentation 

* Steady-state 
Probe 

* Dynamic Probe 


(a) Inlet Test Model: 


(b) Total Pressure: 


Ap t (i,j) fluctuating 



steady state 


Time-sec 


p t (i ’ j) - p s £ 

Nearly Normal 
Distribution 


(c) Distortion Factor (resulting from a combination of all probes) 


DYNADEC Peak 


Time-sec 


instantaneous peak distortion 
factor. 


'M.nst 
Nearly Normal 
Distribution 
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In the second category, the peak dynamic distortion is 
predicted using some form of statistical analysis (refs. 3, 

4, 5, 6, 7, 8) which is more economic than the DYNADEC 
approach. There are three major statistical approaches. 

The first one, due to Jacocks, et al., uses Gumbel’s extreme 
value theorem to extrapolate the peak dynamic distortion 
within a specified time period from the first short time 
segment of deterministic information (ref. 4). The second 
approach, due to Motycka, utilizes a random number generator 
to simulate the total pressure fluctuations at each probe 
based on the measured steady state total pressures and the 
measured total pressure rms levels. Time variant distortion 
factors are calculated using the generated instantaneous 
total pressures, from which the peak dynamic distortion is 
determined (ref. 5). The third approach, developed by 
Melick, et al., is a modeling approach that postulates that 
the dynamic effects in the inlet flow is totally attributed 
to the pressure disturbances caused by a sequence of 
convecting vortices in the inlet duct. The peak dynamic 
distortion is determined from this physical model and its 
statistical properties. Filtered and Unfiltered total 
pressure rms levels are used to identify main variables of 
the vortex flow model (ref. 3, 6). Among all the methods in 
this category, the Melick method is the most efficient in 
terms of numerical analysis whereas it is not as accurate as 
the Motycka method in predicting the peak dynamic distortion 
patterns (ref. 7). 
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Since the forementioned approaches require measurements 
in the time variant total pressures using high response 
total pressure probes, design of high performance inlet is 
costly in terms of instrumentation, analysis and testing, 
time. Thus, it is highly desirable to develop analytical 
means for predicting peak dynamic distortions without 
dynamic measurements. 

There are analytical methods of predicting steady state 
distortions but not the dynamic peak distortions. Also, 
these analytical methods are limited to qualitative 
predictions only, since the total pressure losses through 
the inlet duct can not be well predicted. A summarized 
review of the current analytical methods is included in 
reference 9. Reference 10 provides an application of the 
time marching finite volume method to the prediction of 
steady, quasi-invicid inlet duct flows. 

The present work makes the first attempt ever to 
predict the peak dynamic distortion based on the steady 
state total pressure measurements at the compressor face 
station. Once the approach of the present work is 
successful, it would be highly encouraging to extend the 
present analysis to a pure dynamic distortion predicting 
technique in the future research, without any experimental 
measurement, by incorporating the present analysis in an 
existing steady state distortion prediction code. This 
would allow early identification of possible dynamic 
distortion problems before the first wind tunnel model is 
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built. 


In the present work, an analytical method of predicting, 
the dynamic peak distortion, based on the steady state total 
pressure measurements, is proposed and evaluated. First, 
axial velocity distributions of the inlet flow are 
calculated from the measured steady state total pressures. 

A two-equation turbulence model, k-e model, is then employed 
to solve for turbulent kinetic energy distributions over the 
compressor face using the calculated velocity distributions 
and statistically estimated boundary conditions. The total 
pressure rms levels are computed from the predicted 
turbulent kinetic energy. Finally, the dynamic peak 
distortions are estimated by the Melick program using the 
predicted total pressure rms levels. 

In solving the turbulent equations which are coupled 
and nonlinear, a finite difference scheme is used to 
discretize these equations such that tridiagonal matrix 
equations are obtained for successive line relaxation 
iterations. Effects of the relaxation factor on the 
convergence of the numerical scheme are examined from which 
an optimum computational scheme is determined. 

Sensitivity study of the variables and parameters of 
the turbulent equations to the final results is also 
included to reveal how the solutions vary with the variables 
and parameters. 

Results of the present method are compared with those 
of the DYNADEC for five inlet data sets. These data sets 
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cover a range of Mach numbers from 0.6 to 2.5, a range of 
angles of attack from -10 degree to 25 degree and a range of 
angles of sideslip from 8 degree to -10 degree. Inlet 
configurations involved in these investigations are 
illustrated in Figures 1 through 5. 

Major objectives of the present study are: (1) to 
provide an efficient tool for the preliminary inlet design 
and development in which only the steady state total 
pressure measurements are required; (2) to demonstrate that 
the peak dynamic distortions can be predicted by using the 
k-e turbulence model; and (3) to suggest possible future 
studies in improving the analysis of the inlet flow dynamic 
distortion. 
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3. METHOD OF APPROACH 


In the present analysis, a semi-empirical approach in 
predicting dynamic peak distortion of the inlet flow, based 
on compressor face steady-state total pressure measurements, 
is developed under the following major assumptions: (l) 

turbulent characteristics of the inlet flow, at the 
compressor face station, can be represented by an 
incompressible, high Reynolds number, two-equation 
turbulence model -- k-e model (ref. 11); (2) the inlet flow 

near the compressor face is dominated by the axial velocity 
component such that secondary flow (flow components that are 
not parallel to the inlet duct) effect can be neglected; 

(3) streamwise variation of the turbulent characteristics is 
small and can be ignored from the turbulence model; and (4) 
boundary (near-wall) total pressure rms level is 
proportional to the total pressure loss near the wall of the 
inlet duct so that a statistical correlation between the 
measured steady-state total pressures and the boundary total 
pressure rms level can be developed using the HiMAT data 
set . 

Assumption (l) is justifiable based on the following 
reasons: (a) compressibility effect on the turbulent 

characteristics of the inlet flow can be neglected since the 
flow near the compressor face station is always at 
low-subsonic speed; (b) the k-e turbulence model is well 
tested for many types of turbulent flow such as turbulent 
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boundary layer along a flat plate (refs. 11, 12, 13, 14, 15, 
16), turbulent flow inside rectangular channel or circular 
pipe (refs. 17, 18, 19, 20), free-shear turbulent flows 
(refs. 21, 22, 23) and other applications (refs. 24, 25); 

(c) although the accuracy of the k-e model is comparable to 
other two-equation turbulence model, such as the one 
developed by Rotta and Vollmers (refs. 26, 27), it is found 
that the k-e model is more efficient numerically. 

Despite the fact that assumptions (2) and (3) may not 
be realistic from a theoretical point of view, for some 
inlet configurations and flight conditions, these two 
assumptions are imposed because there is normally no flow 
direction or streamwise turbulence information available. 

It will be shown later that assumption (4) makes a 
major contribution to the present analysis. Assumption (4) 
is made in order to develop a reliable method for estimating 
the boundary conditions for the turbulent model. 

Note that the present analysis does not take into 
account the effect of boundary layer bleed immediately 
upstream of the compressor face since turbulent structure of 
the inlet flow would be quite different in that case. 
Therefore, it is expected that the present model may not be 
suitable for inlet flow with strong boundary layer bleed. 

Based on these four assumptions, developments of the 
mathematical formulations and the method of estimating the 
boundary conditions are described in the following two 
sections . 
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3.1 Mathematical Formulations 


Partial differential equations of the k-e turbulent 
model suitable for incompressible high Reynolds number 

turbulent pipe flow (ref. 12) can be written as: 

C C 

U 1 / _U^_ 

sr r 39 rsr^e 3r ; "Fa© v e 


+C u f[(f|) ,+ -F<§§>‘l-* 
c c 

_ 3 e._ 3 e.wae lj , uk* ae n , 1 3 / u^3e u 

u H +v 5F + 7^ 7i? ( dTr r iF )+ F39 ( ^ra9 ) + 


(1) 


( 2 ) 


where u, v and w are time-averaged velocity components in 
axial, x, radial, r, and circumferential, 9, directions 
respectively, k represents turbulent kinetic energy, e is 
turbulent kinetic energy dissipation rate, is radial 

a "C 

velocity gradient, is circumferential velocity gradient 
and the five empirical constants are: 
d k - 1.0 

a £ = 1.3 
C = 0.09 

U 

c el * i' 45 
C c2 = 2,0 

Detailed derivations of equations (l) and (2) are included 
in Appendix A. 

Using assumptions (2) and (3), the left hand side of 
equations (l) and (2) can be eliminated. This results in: 
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1 _3_ / C U k* 3k s 
far v O^e ar ’ 


,1 a 
+ Fae 


/ ii k 2 £k \ 
( ^e 3© 



£ = 0 


(3) 


o c 

1 ) / uk? 3e \ , 1 5/ u 3 G \ . £• r iff (liLy +. 1^ /ju ui ^ s_ n 

fif ( f“r' r ^ )+ Fae ( "ae )+C el C u kl V + F% ;) J C e2k 0 


(4) 


In equations (3) and (4), the radial and 
circumferential velocity gradients are required as inputs 
which can be obtained from the measured steady-state total 
pressures. First, the axial flow velocity at each pressure 
probe location is computed from the steady-state total 
pressure measurement by the following equation: 



where is face-average flow velocity, p t and p g are total 
and static density respectively, and P g are total and 
static pressure respectively and 7 = 1.4. Derivation of 
equation (5) is described in more detail in Appendix B. 

A cubic spline interpolation subroutine is then 
employed to calculate and interpolate the radial and 
circumferential velocity gradients at every computational 
nodal point on the compressor face. An artificial smoothing 
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routine is also used to account for large velocity 
variations near the wall. 

After the solution of turbulent kinetic energy is 
obtained from equations (3) and (4), total pressure rms 
level is computed from the turbulent kinetic energy using 
equation (6). 


rms 


Pt (7 - 1 ) 



J( 2tf k/3) + 


( 6 ) 


Note that rms has been normalized by the face-average total 
pressure (in psi), P^. Appendix C describes the detailed 
derivation of equation (6). This equation is only an 
approximation since the relation of isotropic turbulence has 
been assumed in the derivation. 

3.2 Boundary Conditions 

In this section, a statistical correlation between the 
boundary total pressure rms level and the steady-state total 
pressure is obtained by applying assumption (4) to the HiMAT 
data set. Assumption (4) can also be expressed in the 
following equation: 


°k = 1 - (P 


(rms) 


nre - o - average _ 


nnSr 


t2 


1 - ( P h0 V 

ring - o - average t 2 o 


(7) 
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where rms, = (rms) c , 

b x 'nng-5-average ’ 

^t2^b ~ ^ P t2' > ring-5-average’ _ 

C k stands for a correlation factor and C P 1 2 ^ t> an< ^ rms b ^ ave 
been normalized by § t g and P^ respectively. P t g and P^ 
are free stream and face-average total pressure respectively 
(in psi). 

Obviously, is not a constant due to the fact that 
different turbulent structures of the inlet flow are 


produced in different flight conditions. Thus, a different 
proportion of the turbulent kinetic energy in the total 
energy loss is expected for different turbulent structures. 
Nevertheless, by examining the HiMAT inlet data set, there 
is a fairly consistent relation between the turbulence level 
and the face-average total pressure loss, (l-P^). From 
this clue, a statistical analysis between and the 
face-average total pressure loss is done from which a second 


correlation between and (1-P is extracted from the 


HiMAT data set, 


Figure 6 shows a linear regression of log^gC^ on 
l°glo(l-Pt2^ ^ or t ' ae i n l et data at subsonic and 

transonic speed. Data at supersonic speed will be handled 
differently, as will be discussed later, for the effect of 
inlet shock system. From Figure 6, the equation of a 
regression line is obtained. This is expressed by equation 
( 8 ). 
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log 10 C k = -°* 9393 + (-0*2537) log l0 (l - P t2 ) 


( 8 ) 


or, 


= 0.115 (i - p^r 0,2587 


(9) 


It is shown clearly from Figure 6 that the regression line 
has a negative slope. This reveals the fact that the higher 
the turbulent level of the inlet flow the higher the total 
energy that is dissipated through turbulent mixing and 
redistribution up stream of the compressor face. 

Tests of this regression line on other inlet data is 
also illustrated in Figure 6. It can be seen that this 
regression line, equation (8), represents the inlet flow 
turbulent characteristics of other inlet configurations 
equally well. 

The same analysis is applied to supersonic cases of the 
HiMAT data set without accounting for the effect of inlet 
shock system. This is shown in Figure 7, which reveals 
discrepancy between equation (8) and the test data. The 
reason for this discrepancy is that the total pressure loss 
through the inlet shock system is due to the shock 
compression process which is not responsible for the major 
part of the generation of turbulent kinetic energy. 

In order to obtain a consistent correlation between 
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and the steady-state total pressure for supersonic cases, 
equation (7) is revised to be: 


°k 


tl 


msb 

^T"P 


1 2 'h 


( 10 ) 


where P^, normalized by the free stream total pressure, 
stands for the total pressure recovery through the inlet 
shock system. If P ^ is not measured in the test, it can be 
estimated from a set of curves of maximum total pressure 
recovery for 2-dimensional oblique shock and conical shock 
systems as shown in Figure 8 (ref. 28). 

Again, a correlation between and the face-average 
total pressure loss after inlet shock system, (P^-P^K ^ or 
HiMAT inlet data at supersonic speed, is shown in Figure 9. 
Corrections for the effect of inlet shock system are 
obtained from the curve of one 2-dimensional oblique shock 
in Figure 8. A revised version of equation (8) is also 
shown in Figure 9. Equation (ll) is the revised version of 
equation (8). 


log 10 Ck = -°- 9393 + (-0-2587) log l0 (P tl - P^) 


( 11 ) 
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or, 


- 0-U5 ( p t r P t2> 


-0.2587 


( 12 ) 


In Figure 9, equation ( 1 1 ) represents higher value in 
than the data points. This is due to the fact that the 
shock correction curves given in Figure 8 is the least 
amount of correction that a inlet can possibly attain under 
any flight condition. In reality the total pressure 
recovery through an inlet shock system is lower than that 
obtained from Figure 8 such that the mean of the data points 
in Figure 9 would be very close to equation (11). However, 
it is feasible to use equation (ll) and Figure 8 for the 
estimation of since a slight over estimation of can be 
attained. This means that the boundary rms level would be 
slightly over estimated which is desirable for a 
conservative design process. 

The validity of equation (11) is tested using other 
inlet configuration of Melick's supersonic test cases as 
shown in Figure 9. The correction for the effect of the 
inlet shock system is also accounted for by using Figure 8. 
Good correlation between equation (11) and the test data is 
shown clearly in Figure 9. 

Conclusion of the statistical study described in this 
section can be summarized by the following equations which 
are combinations of equations (7) and (9), and equations 
(10) and (12). 
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For subsonic and transonic cases: 


= 0.115 (1 - P t2 r a2587 . [1 - (13) 


For supersonic cases: 


- 0.115 <P U - P t2 )-°' 2587 . [P tl - (P t2 ) b ] 


(14) 


Also, equation (6) can be written as: 


k 




3 

rms t2 'Y u 


(15> 


Using equations (13), (14) and (15), boundary conditions for 
the turbulent kinetic energy, k, of the k-e model can be 
estimated. 

Note that the estimated boundary rms level is an 
average value near the wall rather than a circumferential 
distribution as the actual measured rms level would be. It 
is not only that there is no consistent rule for estimating 
the circumferential rms level distribution but that the use 
of an average boundary rms level predicts almost the same 
face-average rms level as that predicted by using the exact 
boundary rms level distribution. This will be described 


16 


later in more detail in section 5.2. 

Boundary condition for the turbulent kinetic energy 
dissipation rate, e, is also required in order to solve 
equations (3) and (4). This boundary condition can be 
obtained by applying a wall function approach (ref. 11, 12, 
17, 18) to the k-equation, equation (3). The wall function 
approach states that the turbulent energy generation term, 
the third term of equation (3), and the turbulent energy 
dissipation term, the last term of equation (3), are almost 
equal near the wall. The following equation supplements 
this statement. 


e k ] ( -W i+ T T< %5 )1 (16) 

Using equation (16), the boundary condition for the 
turbulent kinetic energy dissipation rate can be estimated. 

3.3 Interfacing with the Melick Method for Predicting 
Peak Dynamic Distortions 
Besides the compressor face steady-state total 
pressures and the face-average total pressure rms level, a 
ratio of the filtered to unfiltered mean square (ms) levels 
and the corresponding cut-off frequency are required as 
inputs for the prediction of dynamic peak distortion factors 
by using the Melick method (refs. 3, 6). Unless this ms 
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ratio and the cut-off frequency are determined 
experimentally, an approximation is suggested here. It is 
learned by examining, the available test data sets that the 
ratio of filtered to unfiltered ms levels is around 0.5 at a 
cut-off frequency ranging from 800 hz to 1000 hz . This 
corresponds to fairly constant frequency response of the 
turbulent characteristics of the inlet flow. Also, it is 
mentioned in reference 3 that the prediction of dynamic peak 
distortion is not very sensitive to the ms ratio and the 
cut-off frequency. In the present analysis, an ms ratio of 
0.5 at a cut-off frequency of 1000 hz are used for every 
data set. It will be shown from the results that this 
approximation is reasonable. 
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4. NUMERICAL SCHEME 


Since equations (3) and (4) are coupled and nonlinear, 
a finite difference numerical scheme is used to discretize 
these two equations (refs. 29, 30, 31). A computational 
grid system and a set of finite difference equations are 
arranged in such a way that an implicit numerical scheme is 
obtained. A set of tridiagonal matrix equations are also 
attained for efficient successive iterations. A relaxation 
factor and grid size parameter are determined in the 
analysis of numerical convergence in section 4.3. 


4.1 Grid System 

A grid system based on the locations of total pressure 
probes on the compressor face is illustrated in Figure 
10(a). Figure 10(a) shows a 40-probe compressor face 
instrumentation configuration. For better finite difference 
approximation, finer grids are used by dividing the space 
between probes into N equal spaces. N , an integer, is 

§ o 

designated as the grid size parameter. In this way, smaller 
grid size can be attained near the wall where larger 
variation in turbulent characteristics is expected (ref. 

18). Total number of non-boundary grid points of the grid 
system can be calculated by the following relation. 


TO = (N Y. N. • N p 


(17) 
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where TN denotes the total number of non~boundary grid 
points, N is the grid size parameter, N and N represent 

o ^ r 

the number of rings and the number of rakes of the total 
pressure probes respectively. 

4.2 Finite Difference Formulations 

In this section, a set of discretized difference 
equations are derived by applying central difference scheme 
to equations (3) and (4). Referring to Figure 10(b), the 
central difference equations, in radial (r) and 
circumferential (0) directions, for any variable 4> about 
point p can be written as (ref. 30): 


3$ 

3r 


Ar. 


Ar^i+ATi 


£ 


Ar i-1 


Id) 


^i- 


Ar i- 


ATi 


(_!i±Lil^l) (18) 



2 A6 


(19) 


Note that an irregular grid spacing has been assumed in the 
radial, r, direction. 

Applying equations (18) and (19) to equations (3) and 
(4), by changing $ to k and e respectively, the following 
finite difference equations can be obtained. 
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A- i k. . .+ A. ~ k. .+ A. , k.^, • = C. 
"x,l l-l, j 1,2 i,j 1,3 i+l, j l 


B. 1 e. 1 .+ B. t e. .+ B. , e. ... . = D. 
1,1 l-l, j i,2 i,j 1,3 i+l, j l 


(21) 
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Note that equations (3) and (4) can not be applied to 
the point at the center of the inlet duct where r=0. This 
problem can be handled by rewriting equations (3) and (4) in 
Cartesian coordinates. The central differencing technique 
is then applied to the new equations about point p and the 
surrounding points n, e, s and w shown in Figure 10(c). 

Equations (20) and (21) are then written for all points 
on a radial line, as shown in Figure 10(a), from point 0 to 
point n+1. From this, two tridiagonal matrix equations are 
obtained for efficient line relaxation iterations. These 
two matrix equations are expressed as: 
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Equations (22) and (23) can be solved easily by the routine 
described in Appendix D. 

In order to start the implicit iterative numerical 
scheme, initial guess of the values of k and e for every 
non-boundary grid point must be given. Initially, the 
boundary values of k and e, estimated from equations (13), 
(14), (15) and (16), are assigned to every grid point to 
start the iteration. 

A relaxation factor, f r , is then used to control the 
convergence of the numerical scheme (ref. 3l) • Effect of 
the relaxation factor on the rate of convergence will be 
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presented in the next section. Using the relaxation factor, 
f r , the updated solutions of k and e are obtained from the 
following relations: 

k (wl> _ k (n) + £__ [ k’ - k (n) ] (24) 

c (n+l) - c (n) + £ r [ c* - e (n) ] (25) 

where k" and e" are solutions of equations (22) and (23), 
and superscripts (n) and (n+l) stand for the old and new 
values respectively. 

4.3 Numerical Convergence 

Numerical convergence of the finite difference scheme 
described in the previous section depends on two factors: 
the relaxation factor, f r , and the grid size parameter, N . 
To study the effect of f r and N^ on the rate of convergence 
of the numerical scheme such that an optimum numerical 
scheme can be determined, a relative percent error of each 
iteration is defined as: 


E 


r 



+ 




] • 100 % 


(26) 


where TN is the total number of non-boundary grid points and 
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the sign, ]T, stands for the summation over all non-boundary 
grid points. 

Steady-state total pressure data used in this study is 
obtained from the HiMAT data set. 

To determine the value of f r , the relative percent 
error, E r , at the end of the 30th iteration for several 
values of f and boundary rms levels, rms-^, is computed and 
presented in Figure 11. It is seen clearly from Figure 11 
that E r is decreasing with the increase in f , until 
exceeds a boundary beyond which the numerical scheme becomes 
divergent. The boundary for divergent and the rate of 
convergence depend also on the boundary rms level, rms^. 

From this a value of f r =1.0 is selected for relatively good 
characteristics of convergence and allowing a margin to the 
divergent boundary for safety operation of the numerical 
scheme . 


To determine the grid size parameter, N , the relative 

§ 

percent error in 30 iterations, E^, for several values of 
Ng is computed and illustrated in Figure 12. f r and rms^ 
are fixed at values of 1.0 and 0.011 respectively. It is 
shown in Figure 12 that E stays at almost a minimum value 
beyond N^=4. For efficient numerical scheme in terms of 
computing time and good characteristics of convergence, N =4 

5 

is chosen for successive computations. 
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5. SENSITIVITY STUDY 


Prediction of the face-average total pressure rms level 
and the dynamic peak distortion by using the present 
analysis depends on the following factors: (l) the measured 

steady-state total pressure profiles; (2) the estimated 
boundary rms level; and (3) the empirical parameters of the 
k-e turbulence model. In the following sections, effects of 
these factors on the prediction will be examined from which 
the sensitivity of these factors on the present approach can 
be identified. 

5.1 On the Velocity Distributions 

Study included in this section will concentrate only on 
the effects of radial variations of the flow velocity. 

In Figure 13(a), five radial steady-state total 
pressure profiles, which are equivalent to the velocity 
profiles, are assumed. These total pressure profiles cover 
a range of boundary layer characteristics from a favorable 
boundary layer, profiled); to a separated boundary layer, 
profiled). Note that prof ile (4) stands for the boundary 
between separated and non-separated boundary layer profiles. 
To predict the dynamic distortion, in this study, a constant 
boundary rms level, rms^=0.02, is used for each pressure 
profile. 

Results of the prediction using the present analysis, 
based on the assumed five steady-state total pressure 
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profiles are illustrated in Figure 13(b). In Figure 13(b), 
only prof ile (I) produces lower rms level near the center of 
the inlet duct. The rms level is heightened near the duct 
center as the total pressure profile is close to the 
separated flow. 

Figure 13(c) shows the predicted face-average rms 
level, rmSp, and the dynamic contribution to a distortion 
factor, K^ 2 > for the five total pressure profiles. It is 
shown clearly that the profile with a favorable boundary 
layer, profile (T), produces the lowest rms^ and dynamic 
contribution to while the separated flow, profile(5X 
predicts the highest rms^ and dynamic contribution to K A2- 

This study shows that: (l) the present method is 

sensitive to the measured steady-state total pressure 
profiles; and (2) the dynamic contribution to the peak 
distortion can be minimized by designing the inlet contour 
or using boundary layer treatment devices, (e.g. vortex 
generators,) so that a favorable boundary layer profile can 
be attained. 


5.2 On the Boundary Conditions 

In this section, effects of the distribution and 
magnitude of the boundary rms level, rms^, on the prediction 
are presented. Steady-state total pressure data used in 
this investigation is obtained from an inlet test case, with 
the test condition shown on the top of Figure 14 (data point 
1948,) of the HiMAT data set with vortex generators. 
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First, the effect of the boundary rms distribution on 
the predicted rms distributions on the compressor face is 
illustrated in Figure 14. In Figure 14, the circular symbol 
stands for the measured rms levels, the solid line is the 
result of the prediction using a constant boundary rms level 
estimated by equation (13) and the dashed line represents 
the result of the prediction based on the actual measured 
boundary rms level, which is not a constant 
circumferentially. The indicated number of degrees on the 
upper left corner of each of the eight figures in Figure 14 
denotes the angular position of the pressure rakes on the 
compressor face. The measured and predicted face-average 
rms levels, rms, are also shown in the upper portion of 
Figure 14. 

It is shown clearly in Figure 14 that the use of the 
actual rms-^ distribution improves the prediction, especially 
for the pressure rake at 180 degree. However, a discrepancy 
in the prediction of face-average rms level, rms , is not 
very significant. 

In Figure 14, discrepancies between the measured and 
predicted rms distributions can be seen clearly. This is 
due mainly to the inevitable assumptions of the theoretical 
model, assumptions (2) and (3), in which the effects of the 
secondary flow and the axial variations of k and e are 
omitted. 

Note that the computed value of rms using the constant 
rms-^ is slightly over-predicted over the measured rms and 
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the one computed using the actual rms^ distribution is 
under-predicted. The same result is obtained for several 
test cases of the HiMAT data set. From this, the use of the 
estimated constant rms-^ in the present analysis is feasible 
in the inlet design process. 

Next, effects of the magnitude of rms^ on the predicted 
value of face-average rms level, rrhs ^ , and the dynamic peak 
distortion factor, p ea ^> are presented in Figure 15(a) 

and Figure 15(b) respectively. It is shown in Figure 15 
that rmSp and p ea k are ver y sensitive to the variation 

of rms^. This means that good predictions of the 
face-average rms level and the dynamic peak distortion 
depend on a good estimation of the boundary rms level. 
Validity of the present analysis may largely be attributed 
to the statistical estimation of the boundary rms level, 
rms-^. 


5.3 On the Parameters of the Turbulence Model 

Effects of the values of the five empirical constants, 
0^, d £ , 0^, and C^, on the predictions of face-average 

rms level, rms^, and dynamic peak distortion factor, 

p ea k> are investigated in this section. This is done by 
varying one of the parameters while keeping the other 
parameters at constant values. Results of this study is 
illustrated in Figure 15. 

It can be seen clearly from Figure 15 that the 
predicted face-average rms level and the dynamic peak 
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distortion factor are not sensitive to the changes in the 
five parameters of the turbulence model. Therefore, an 
accurate estimation of the boundary rms level is far more 
crucial than any adjustment in the empirical constants of 
the turbulence model. 
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6. NUMERICAL RESULTS AND DISCUSSIONS 


In the following sections, numerical predictions by the 
present analysis and data comparisons with five inlet data: 
sets are included. The five inlet data sets consist of 
experimental results of subsonic, transonic and supersonic 
inlet configurations under various flight conditions. Inlet 
configurations and some measured results of these data sets- 
are illustrated in Figure 1 through Figure 5. Three aspects 
of data comparisons will be included. They are; 
face-average total pressure rms level comparison; peak 
dynamic distortion factors comparison; and detailed rms 
level distributions on the compressor face comparison-. 

6.1 The HiMAT Inlet Data Set 

Reference 34 describes an experimental investigation of 
a subscale HiMAT model with forebody, canard, inlet duct and 
optional boundary layer control devices, vortex generators. 
The HiMAT inlet tests were conducted by NASA in the NASA 
Lewis 8'x6’ supersonic wind tunnel facility. The HiMAT 
model has an under-fuselage inlet with a high-divergence 
S-shape subsonic diffuser. At the compressor face station, 
there were 40 steady-state total pressure probes and 40 
dynamic high response total pressure probes installed in an 
8-rake, 5-ring configuration as illustrated in Figure 1. 

Using the present analysis, face-average total pressure 
rms levels and peak dynamic distortions are predicted based 
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on the steady-state total pressure measurements of the HiMAT 
data set. Comparisons of the analytically predicted and 
experimentally measured rms levels and peak dynamic 
distortions, for the HiMAT inlet model without vortex 
generators, are shown in Figure 16. Reasonably good 
accuracy of the present method in predicting the 
face-average rms levels is revealed from Figure 16. 

However, Figure 16 also shows that the present analysis 
underpredicts the peak dynamic distortion factors for most 
cases of the HiMAT data set without vortex generators. 

Since these test cases contain incipient or intermittent 
separated flow, under-prediction of the Melick statistical 
method can be expected (refs. 3, 34). Note, in Figure 16, 
the measured peak dynamic distortion factors, APRS m (defined 
in Table 1,) are the DYNADEC results (ref. 1). 

Figure 17 illustrates comparisons of the peak dynamic 
distortions, in terms of IDC , predicted by the present 
analysis and the Melick method based on measured rms levels 
for the HiMAT data set without vortex generators. It is 
shown in Figure 17 that predictions of the present analysis 
compare well with the Melick predictions. Note that the 
peak dynamic distortion factors of Melick' s prediction, 


IDC 


max, Melick' s-prediction 7 


are obtained from the Melick 


statistical predictions using 40 total pressure rms 


measurements . 


Figure 18, (a) through (c), presents comparisons of the 
predicted and measured rms level distributions at the 
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compressor face for three test cases selected from the HiMAT 
data set (without vortex generators.) Discrepancies between 
the predicted and measured rms level distributions can be 
observed in Figure 18. Since the HiMAT inlet model has a 
short S-shape inlet duct with large streamwise curvature, 
substantial secondary flov; components can be expected (ref. 
17). The presence of the secondary flow violates 
assumption (2) of the theoretical model. Therefore, some 
discrepencies in the prediction of the rms level 
distributions should be expected. However, the face-average 
rms level is predicted well which is responsible for the 
good prediction of peak dynamic distortion using the Melick 
statistical method. 

In Figure 19, comparisons of the predicted and measured 
face-average rms levels and peak dynamic distortion factor 
are presented for test cases of the HiMAT data set with 
vortex generators. The predictions shown in Figure 19 are 
more accurate than those cases without vortex generators, 
Figure 16. For these cases with vortex generators, boundary 
layer of the inlet flow is controled by the vortex 
generators successfully (ref. 34), which results in better 
predictions of the Melick statistical method. 

Figure 20 illustrates comparisons of the peak dynamic 
distortions, in terms of IDC ma „, predicted by the present 
analysis and the Melick method based on measured rms levels 
for the HiMAT data set with vortex generators. It is shown 
in Figure 20 that predictions of the present analysis 
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compare well with the Melick predictions. Note that the 
peak dynamic distortion factors of Melick's prediction, 

lDC .ax, Melick’ s-predictio n > are obtained from the MeUck: 
statistical predictions using 40 total pressure rms 

measurements. 

Data comparisons of the rms level distributions are 
illustrated in Figure 21, (a) through (c). These cases are 
selected from the HiMAT data set with vortex generators. As 
for the previous cases, discrepancies in these comparisons 
can be observed. Further study is needed to improve the 
prediction of rms level distributions on the compressor 
f ace . 

6.2 The Melick Test Cases 

Configuration of a supersonic mixed compression 
axisymmetric inlet model, used to supply the test cases for 
Melick's computer program, is illustrated in Figure 2 (ref. 
35). This inlet model consists of a center body with fixed 
cone angles. Free stream Mach numbers of the test are 1.6 
and 2.5 with total pressure recovery through the inlet shock 
system assumed to be 0.995 and 0.9 respectively for this 
type of inlet. Three test cases are available for data 
comparisons . 

In Figure 22, comparisons of the analytically predicted 
and experimentally measured face-average rms levels and peak 
dynamic distortion factor are illustrated for the three test 
cases. Good predictions by the present method are shown 
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clearly in Figure 22. Note the measured face-average rms 
level, rms , is an average of 14 total pressure rms 
measurements and the peak dynamic distortion factor of 

MeLick's prediction, K A2 ,Melick' s-prediction - is ° btained 
from the Melick statistical prediction using the 14 total 
pressure rms measurements. For these test cases, there are 
no DYNADEC results available for data comparisons. 

The compressor face rms level distributions are 
predicted and compared with the measured results. This is 
shown in Figure 23, (a) through (c). Again, discrepancies 
in these comparisons can be seen in Figure 23. Although 
this inlet model is axisymmetr ic , substantial secondary flow 
components can still be expected with an inlet angle of 
attack or sideslip angle. 

6.3 Subsonic Inlet 

Configuration of a subsonic full scale short S-shape 
inlet model, which is very much like the HiMAT inlet model, 
is shown in Figure 3. Center line of the engine is tilted 
as shown in Figure 3. Free stream Mach number of the test 
is subsonic. There are six test cases available for data 
comparisons. These data were provided by the Air Force 
(AFFDL), Wright-Patterson Air Force Base, Ohio. 

Comparisons of the predicted and measured results are 
illustrated in Figure 24, Figure 25 and Figure 26, (a) 
through (f). Reasonably good accuracy of the present 
analysis in predicting the face-average total pressure rms 
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levels and the peak dynamic distortion factors is shown in 
Figure 24. Note the measured peak dynamic distortion 

factor, IDC , is obtained from the DYNADEC system (ref. 

7 max 7 1 v 

m 

1 ). 

It is seen from Figure 24 that the present method 
slightly underpredicts the peak dynamic distortion factor 
for some cases. This is due to the limitation of the Melick 
statistical method that it can not properly predict the peak 
dynamic distortion for inlet flow with separated boundary 
layer (ref. 3). Unfortunately, every case of this subsonic 
inlet test data set contains separated boundary layer at the 
compressor face station. Further study is required to 
improve the Melick statistical method to handle the 
separated flow conditions. 

Good comparisons of the predictions of the present 
analysis and the Melick method based on rms measurements are 
illustrated in Figure 25. 

Similar comparisons of the rms level distributions of 
the six test cases are shown in Figure 26, (a) through (f). 

6.4 Transonic Inlet 

Configuration of a transonic subscale long S-shape 
inlet model is illustrated in Figure 4. Six test cases of 
this inlet model are available for data comparisons. These 
data were also provided by the Air Force (AFFDL) . 

Comparisons of the predicted and measured results are 
shown in Figure 27, Figure 28 and Figure 29, (a) through 
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(f). Good accuracy of the present method in predicting the 

face-average total pressure rms levels and the peak dynamic 

distortions is shown clearly in Figure 27. In contrast to 

the results described in the previous section, or Figure 24, 

the present analysis overpredicts the peak dynamic 

distortions slightly, which is illustrated in Figure 27. 

Since every case of this data set retains attached boundary 

layer at the compressor face station, better prediction of 

the peak dynamic distortions, by the Melick statistical 

method, can be expected. This results in Figure 27. Note 

the measured peak dynamic distortion factor, K A2 > is 

m 

obtained from the DYNADEC system. 

Figure 28 illustrates good comparisons of the 
predictions of the present analysis and the results of the 
Melick predictions based on total presure rms measurements. 

Also, slightly improved predictions of the rms level 
distributions are obtained since a long S-shape inlet duct 
produces less streamwise curvature effect which incurs 
smaller amount of secondary flow components which is closer 
to the theoretical assumptions of the present analysis. 
Results of the predictions of rms level distributions for 
the six transonic cases are illustrated in Figure 29, (a) 


through (f) 


6.5 Supersonic Inlet 

Configurations of four supersonic inlet testing models 
are illustrated in Figure 5. These inlet models are 
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designated as A-l, A-2, B-3 and B-4. There are thirteen 
test cases of the supersonic inlet models available for data 
comparisons. These data were also provided by the Air Force 
(AFFDL) (ref. 36). Data included in this data set are the 


steady-state total pressures and the peak dynamic 


distortions which were screened based on K, 


There is no 


information about the measured rms level available for data 


comparisons. Inlet configurations of the four models are 
described in more detail in the following paragraphs. 

The inlet model, A-l, is a two-dimensional fuselage 
side-mounted external compression inlet. The compression 
surface consists of three variable ramps including, an 
articulated first ramp to allow a variable capture area. 
There is also a variable fourth ramp. A large slot at the 
throat between the trailing edge of the third ramp and 
leading edge of the fourth ramp provides for bypass flow and 
choking control. 

The A-2 inlet design is a 180 degree axisymmetric, 
fuselage side-mounted, external compression inlet. The 
center body consists of a first cone fixed at 18 degree with 
a variable radius second cone that has a range from 12 


degree to 30 degree, 

The B-3 inlet is identical with A-2 except that it is 
located beneath the wing of a flat-bottomed blended body 
fuselage. The subsonic diffuser length is constrained by 
the location of the inlet and engine beneath the wing with 
the resultant decrease in diffuser length and increase in 
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diffuser divergence angle. 


The 8-4 inlet is identical with A-l except that it is 

also located beneath the wing in much the same position as 

B-3 inlet. The diffuser length is constrained by the inlet 

and engine locations as well. 

Results of the prediction of the peak dynamic 

distortion are computed by the present method and compared 

with the measured results. This is shown in Figure 30 for 

IDC and K. 0 . Since the measured peak dynamic distortions 
max i\c- 

of these test cases were screened based on K^j better 
comparisons are shown in Figure 30 for peak K.^* Note that 
most of the test cases are underpredicted slightly. This is 
because most of these test cases have separated boundary 
layers at the compressor face station. This causes the 
Melick statistical method to underpredict the peak dynamic 


distortions . 

Finally, accuracy of the present analysis in predicting 
the face-average rms level and peak dynamic distortion 
factor is summarized in Figure 31. Figure 31 is presented 
in terms of percent error of the prediction compared to the 
measured data of all available test cases of the five 
forementioned inlet test data sets. 

In Figure 31, three histograms of the percent error in 
predicting the face-average rms level and peak dynamic 
distortion factor are presented. If the normal probability 
distribution is assumed for these three histograms then the 
mean and standard deviation for the histogram of rms level 
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are 15.5% and 15.2% respectively. The mean and standard 
deviation for the histogram of peak dynamic distortion 
factor, for test cases without separated flow, are 19.51% 
and 15.57% respectively. The mean and standard deviation 
for the histogram of peak dynamic distortion factor, for 
test cases with separated flow, are -3.31 % and 17.86% 
respectively. The distribution of the rms histogram 
represents the statistical characteristics of the boundary 
condition estimation technique of the present analysis, 
equations (13) and (14). The distributions of the peak 
distortion factor histograms represent the combined 
statistical characteristics of the rms level prediction 
technique and the Melick statistical method. In Figure 31, 
the peak distortion factor histograms show that the present 
analysis has nearly 93% of chance to overpredict the peak 
distortion factor, by about 19.51%, for inlet flow without 
boundary layer separation which is the feature of the 
conventional type of inlet. This is a good design 
characteristics. Nevertheless, for inlet flow with 
separated boundary layer, the present analysis has more than 
55% of chance to underpredict the peak dynamic distortion 
factor, by about -3.31%, which is not a feasible design 
characteristics. Further study is required to improve the 
Melick statistical method to handle the separated flow 
condition . 

This result can be used as an inlet design reference in 
using the present method. In other words, this information 
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can be used in inlet design process in using the present 
analysis to see how accurate the present analysis can be in 
predicting the face-average rms level and peak dynamic 
distortion factors. 

One final comment about the prediction of peak dynamic 
distortion of the inlet flow using the present analysis is 
that the present method, owing to the characteristics of the 
Melick statistical method, will underpredict the peak 
dynamic distortion for separated flow. If it is found from 
the steady state total pressure measurements that the flow 
is separated then the peak distortion factor histogram for 
separated flow, in Figure 31, can be used as design 
reference for the predictions of the peak dynamic distortion 
factors . 
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7. CONCLUSIONS 


An analytical method in predicting total pressure rms 
Level and peak dynamic distortion, based on steady-state 
total pressure measurements, has been developed. The total 
pressure rms level is predicted by a two-equation turbulence 
model -- k-c model. The Melick statistical method is then- 
used to estimate the peak dynamic distortion based on the 
predicted total pressure rms level and the steady-state 
total pressure measurements. 

A finite difference scheme has been used to solve the 
two-equation turbulence model. A statistical correlation 
between the measured steady-state total pressures and total 
pressure rms level near the wall has been developed using 
the HiMAT inlet data set. Thus, boundary conditions of the 
turbulence model can be estimated from the steady-state 
total pressure measurements. 

A central differencing scheme has been applied to the 
turbulence model to discretize the equations of the 
turbulence model. Thus, an implicit line relaxation 
formulation is obtained. The finite difference equations 
have been arranged into a set of tridiagonal matrix 
equations for efficient numerical iterations. 

Numerical convergence of the finite difference scheme 
has been investigated in order to determine a relaxation 
factor and grid size for the numerical scheme. Sensitivity 
study has also been included to see how sensitive the 
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present analysis is to the input data and empirical 
parameters of the turbulence model. 

Accuracy of the present analysis has been demonstrated 
through data comparisons for subsonic, transonic and 
supersonic inlet models at various angles of attack and 
sideslip. Five sets of inlet test data have been used for 
this purpose. Results of the data comparisons have shown 
that the present method predicts good results of the 
face-average rms level and peak dynamic distortion factor. 
Around 20% of accuracy in predicting peak distortion factors 
is attained for subsonic, transonic and supersonic inlet 
test cases without separated flow. Further research is 
required to improve the present analysis to handle the 
separated flow conditions. 

As far as the prediction of peak dynamic distortion 
factor is concerned, the present analysis can be used as an 
efficient preliminary inlet design tool without total 
pressure rms measurements. Only the steady-state total 
pressure measurements are required to predict the peak 
dynamic distortion factor. 

A logical extension of the current work would be to 
predict the dynamic distortion of the inlet flow directly 
from an analytical point of view so that no experimental 
information would be required in the prediction. This would 
require a flow field calculation external to the inlet duct 
and a compressible turbulent flow computation along the 
subsonic inlet duct. This approach would require a large 
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computational effort with long computing time, but it would 
be an inlet design tool that would provide early estimates 
of peak dynamic distortion prior to the time when inlet 
models and test data are available. 
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APPENDIX A 


Derivation of k-e turbulence Model 

The continuity and momentum equations for 
incompressible flow can be written in tensor notation as 
following : 






3P 



+ u 



(Al) 

(A2) 


where u^ and P represent velocity components and pressure of 
the flow respectively. p and u denote flow density and flow 
viscosity repectively. 

For turbulent flow, the velocity components and the 
pressure can be decomposed into a time-averaged component 
and a time dependent component. This is given in the 
following relations. 


u i = u^ + ut^ , P = P + P' 

Substituting these relations into equations (Al) and (A2) 
and taking time averaging of the resultant expressions, the 
following time-averaged equations of equations (Al) and (A2) 
can be obtained. 



0 


(A3) 
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.auf ,p 3 l u, 

P [ 57 = + 3®- (u.U. + u.'u'.)] - - — + u -T-i 
at 3 Xj v ij i u j' J 3x i M 3X‘. 


( A4) 


Subtracting eauation (A3) from equation (Al) and subtracting, 
equation (A4) from equation (A2) , the following time 
dependent equations are obtained. 



(A5) 


AU. 3UI 3U. 3ul 3p , 3*ui 

tr-r+u- — -+u! — - + u!-—^ - - -^-(u. , u'.)] « --£- + u — i 
3t jsx. J3x. j3x ; Jx? 1 y J 3x, M 3xr. 


3X. X 
J 


oX- 
1 


5r j 


(A6) 


Multiplying equation (A6) by uj. and then taking time 
averaging, equation (A7) is obtained. 



(A7) 


Taking summation of equation (A7) for all i=»l, 2, 3 and 
using equation (A5) and the definition of turbulent kinetic 
energy, k-=( l/2)£ uT 1 , equation (A8) is obtained. 


+ u .il£— -2-S& - u’ £-+ k')l + (- TO) S- - 

3t J3X. 3X.-3X- Y p 4 J •** V 3X.' 

J v J — -v - j „ jLi 


"Y 

(I) Diffusion 


(II) Production Viid 


(A8) 


Dissipation 

rata 
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where V = — 
P 


and u! 


u ! 1 


<-r> - «s$* 


have been used. 


The following empirical approximations were suggested by 
Launder (ref. 16). 


(X) sa. -5— i^_) 


“j 


°k sx j' 


«- 4 u j> 


with ■l.O’ 


(III) H C 


where L'. = C^Jk L is the well known Kolmogorov-Prandtl 
expression for turbulent viscosity. Also, from dimensional 
analysis, c can be written as: 



Thus, the following relation is obtained. 



where 0^=0. 09 and Cq= 0.08 are empirical constants. 

Using these approximations, equation (A8) is written 

as : 


Dk 

Et 


_L_ (Si + vifc - (— ) 2 

3x. ^ C ax/ U «. oX. 

J * J J 



(A9) 


It is sheerly empirical that an equation for e was 
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suggested by Launder. It is assumed that the e-equation 
have the same form as the k-equation with different 
empirical constants. This results in the following 
relation. 


, _3_ /_JU kj_i£v f k/ in 2 _ r e-i 

I5E 3x, o e 3x . ' e| - C el Sj e ^3x.' C e2 


9Q. 

2 


(A10) 


J «• 


where d £ =1.3, C £ ^=1.45 and C & 2 = 2-0 are additional empirical 
constants. 

Finally, equations (A9) and (A10) can be rewritten in 
cylindrical coordinates for steady pipe flow: 

Q lk + ^£k + w3k + J. JL ( ) + 

U 3x sr ri? r 3r ^ 0 ^ e sr ' r 30 ^ 


v 


c c 

n ae , -ae , wae JL / V£_ ae \ ,_lJL / Ja£_3£ ) + 

u 5^ + v 5F + 7^ 73? c o~r r iF ;+ F39^o e se ; + 

& ^ 

+ C el C u k ^l? )1+ F ( 5§ )1] - C e2T 


(All) 


(A12) 
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APPENDIX B- 


Measured Steady State Total Pressure 
and Velocity Distributions 

Relationship between the measured steady state totaT 
pressure and the axial mean flow velocity is derived under 
the following assumptions: 

1. The flow is ideal gas. 

2. The flow is adiabatic.. 

3. The flow is steady. 

4. There is no static pressure variation in the 

direction normal to the wall. 

5. Loss in the total pressure measurement is 
accounted for by a recovery factor, r. 

Under these assumptions, total energy of the flow 
without turbulence at the instrumentation station is picked 
up by the total pressure probe. Effect of turbulence will 
be considered later. The following equation describes this 
relation . 

E t = E 1 (Bl) 

where E t stands for the total energy of the flow after the 
flow stagnates and E^ denotes the total energy of the flow 
just before entering the total pressure probe. For 
adiabatic flow, 
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E t = C p T t ’ E 1 = C p T s + 772 + §< curl v > 


where is the specific heat of the flow at constant 
pressure, T fc and T s are total and static temperatures 
respectively and g function stands for rotational energy of 
the flow inside the boundary layer. Assume that g is very 
small and can be neglected since the measurement is not made 
very close to the wall. Also, V denotes the absolute 
velocity of the flow at the total pressure probe location. 
The equations of state of the ideal gas are: 






and 





where '7=1.4. After substituting these relations into 
equation (Bl), the following equation is obtained. 

P P 

-7~T 7T ~-7?r TT" + * V ' t V 1 - uW-h/ (B2) 

T p t ' “s 

Considering loss in the total pressure measurement, a 
recovery factor, r, is introduced for the total pressure 
probe. That is, 

p . 

r » , where P . is the indicated total 

P t 

pressure . 

Equation (B2) becomes: 
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'TT 'ZT 25 7T + % (uI+ v * + 
y-l rp t /-I p s 


( B3 ) 


The effect of turbulence is accounted for by the following, 
relations : 

v v ? i . v v p ; > p t - h + pi • p s * h* pi 


u =■ u + u' , v =■ v + v' , w =* w + w' 

where the pressure, velocity and density terms are replaced 
by their time-averaged terms plus time-variant fluctuation 
terms. Substituting these relations into equation (B2), a 
time-averaged equation and a time-variant equation can be 
obtained. They are: 

h (u 1 * v* + w 1 ) +% (B4) 

p* p' 

^rr^-^TTr + Cau'+ w+ a*') + % (u ,l + v'*+ w’ 1 ) (bs) 

7 p t 1 

Notice that it has been assumed that p' and p’ are small and 
can be neglected in obtaining equation (B5). Introducing 
the definition of turbulent kinetic energy, 
k=(u ,x +v 11 +w 11 ) /2 , and considering the loss in the total 
pressure measurement, equation (B4) results in: 
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Since u>>v and u>>w, the above equation can be further 
simplified by neglecting v and w. 


y P ti 


^ _ jy _ ^s + l-i 

^-1 rp t ~ 7-1 p s * U 


+ k 


or, 



where U 2 is 
Since k/U^ 
by: 


1 


s 2 


■ 7 w ti 


‘ V - 2 


k 

Ul 


(B6 ) 


the compressor face average flow velocity, 
is very small, equation (B6) can be approximated 






V 


(B7) 


where r=1.0 for most of the total pressure probe used in 
experiment. 
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APPENDIX C 


Total Pressure rms Level 
and Turbulent Kinetic Energy 

Relationship between total pressure rms level and 
turbulent kinetic energy is derived in this section by the 
use of equation (B5) derived in APPENDIX B. Assumptions 
stated in APPENDIX B are also employed here. Equation (B5) 
is: 

pi pi 

+ (UU'+VV'-K ww') (u’ l + v'*+ V- 1 ) 

It is further assumed that the fluctuations of the static 
pressure makes only a small contribution in the above 
equation and can be neglected. Taking root mean square 
(rms) on both sides of the above equation, the total 
pressure rms level can be written a3 : 

rms => J?^ z J [(Qu'+ 7v'+ Sw') + h (u' 1 * v'* + w' 1 )]* 

Since u>>7 and Q>>w and neglecting the higher order terms, the 
above expression can be simplified as: 


P t Ol) 

rms ^ — - 

'y 


Jtfu ' 1 + k 1 


(Cl) 


For isotropic turbulence (i.e. for flow away from the wall), 



u' 2 can be written in terms of k (ref. 12). That is: 


u Tr =* (2/3) k 

Substituting this relation into equation (Cl), an expression- 
in relating total pressure rms level, rms , and turbulent 
kinetic energy, k, can be written as: 

p Cr-1) f 

J( 2/3)u l k + k 2 (C2> 

Equation (C2) is an approximation for flow which is not very 
close to the wall. 
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APPENDIX D 


Solution of Tridiagonal Matrix Equation 


A typical tridiagonal matrix equation can be writ 


as : 


AX =■ F 


where 


a l C 1 ^ 


\ 


\ 


\ 

b 2 a 2 c 2 N 


\ 


\ 


\ 

* \ 


V 


\ 


\ 


\ 


\ * 


\ b 
\ 


\ 

• \ 

n-1 a n-l c n-l 


\ 


\ 


\ b 
\ n 


n 


f x ' 
X 1 


X » \ 


F » < 




The tridiagonal matrix, A, can be decomposed into two 
matrices, (J and V. That is: 


A = UV 


ten 


(Dl) 


(D2) 
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or, 


*1 °l\ 

\ 0 
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\b 3 a 3 c.^ 

> ^ • • N 
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\1 


where p 1 » a 1 , = c 1 /p 1 

p i 3 a i~ b i r i-l’ 1 = 2 » 3 » 4 > » n 

r i s C i/Pi> i = 2, 3, 4, , n 

Using relation (D2), equation (Dl) can be written as: 

UVX = F 

Letting VX ** G, the above relation becomes: 

UG = F ( D 3 ) 

Thus, solution for G can be obtained from the following 
recursion formula: 

Si - V P X 

§i = (f i - b iS i . 1 )/Pi> i = 2, 3, 4, n 
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|)cf In 1 Cion of distortion factors 



























(a) HiMAT Inlet Test Model: 



180° 

(b) Range of Test Conditions: 


Mach No. ** 0.4 — 1.36 
Angles of Attack = -10° - — * +25° 
Angles of Sideslip = 0° — ' -10° 
With and Without Vortex Generators 


Figure 1. Illustration of HiMAT Inlet Test Model and Test Conditions 
(ref. 34) 







(a)Supersonic Mixed Compression Axisyrrmetric Inlet Model: 




(b) Test Conditions and Measured rms of 3 Melick's Test Cases: 


Case No. 

Mach No. 

ot 

P 

rms 

(14-probe average) 

1 

2.5 

5 

0 

.0126 

2 

1.6 

5 

1 

.0194 

3 

2.5 

5 

0 

.0173 


Figure 2. Illustration of the Inlet Test Model of Melick's Test Cases 
and some Test Results (ref. 35) 
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Data pt. 

Mach No. 

P t2 

rms 

tn 

IDC , 

max, peak 

20.40 

subsonic 

.887 

.0336 

.225 

54.30 

subsonic 

.853 

mrvmm 

.326 

81.40 

subsonic 

.925 

.0337 

.127 


subsonic 

3S8 - 

.0537 

.319 



mmum 

.0475 

.329 


subsonic 


.0360 

.144 


Figure 3. Illustration of a Subsonic Inlet Test Model and some Test 
Results (unpublished data from Air Force Flight Dynamics 
Laboratory, Wright-Patterson Air Force Base, Dayton, Ohio) 
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(a) Transonic .15 scale Inlet Model: 



(b) Test Conditions and Seme Measured Results: 


Data pt. 

Mach No. 

P t2 

m 

^A2, peak 

464.12 

transonic 


WKssaam 

IHRKSBHHi 

■ESE1W 

transonic 

mmm 

.0281 


WOEB&m 

.Mm 

mmm 

.0217 H 



transonic 

IliW 

.0414 

1 1 

■MB 

transonic 

wm-Mtm 



wmmvm 

transonic 

BDM 


■BfU MH 


Figure 4. Illustration of a Transonic Inlet Test Model and seme Test 
Results (unpublished data from Air Force Flight Dynamics 
Laboratory, Wright-Patterson Air Force Base, Dayton, Ohio) 
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(a) Configurations of four .25 scale Tailor-Mate Model: 



(b) Test Conditions and Some Measured Results: 



Figure 5. Illustration of four Supersonic Inlet Test Model and some 
Test Results (ref. 36) 
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lliMAT Data Set: Other Available Inlet Data: 

Symbol Mach no. Symbol Mach no. 

a 0.6 A subsonic 


ORIGINAL PAGE IS 
OF POOR QUALfTY 
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Figure 6. Linear Regression of log^C^ on log^Cl-P^) for HiMAT Inlet Data at Subsonic and Transonic 
Speeds 




yrnbol Mach no 


ORIGINAL PAGE IS 
OF POOR QUALITY 



Regression Equation on HiMAT Iniet Test Data pt Superspnic Speed 
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Figure 8. Maximum Total Pressure Recovery of Oblique and Conical Jnlet Shock Systems (Ref. 28 ) 



995 
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Figure 9. Linear Regression of log 1( X, on log 1 r>(P f i -P,.?) for Inlet T§st Data at Supersonic Speed 
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E 

r 



grid size parameter, 

O 


Figure 12. Effect of Grid Size on the Accuracy of Computation 


(a) Assured four total pressure profiles: 



(b) Resultant rms distributions: 



(c) Comparisons of Dynamic Contributions: 


Profile 



rmSp 

Dynamic Contribution 
^A2,peak ^A2,ss^ 

.0195 

.1720 

.0217 

.2382 

.0257 

.3407 

.0283 

.4880 

.0326 

.6688 


Figure L3. Effect of Radial Total Pressure Distributions (equivalent to 
velocity distributions) on the Predicted rms Distributions 
and Dynamic Peak Distortions 
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rms 

m 



Figure 16. Comparisons of the Predicted and Measured rms Levels and 
Peak Distortion Factor for HiMAT Inlet Model without 
Vortex Generators (ref. 34, see Figure 1) 
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Figure 17. Comparisons of peak distortion factors predicted by the present 
analysis and Melick's statistical method based on the total 
pressure nns measurements for HiMAT inlet data set without 
vortex generators (ref. 34, see Figure 1) 
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Figure 20. Comparisons of peak distortion factors predicted by the present 
analysis and the Melick mehtod based on the total pressure rms 
measurements for the HiMAT data set with vortex generators 
(ref. 34, see Figure 1) 
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Figure 22 . Comparisons of the Predicted and Measured nns Levels and 
Peak Distortion Factor for Melick's Test Cases 
(ref. 35, see Figure 2) 
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Figure 24. Comparisons of the Predicted and Measured rms Level and 
Peak Distortion Factor for the Subsonic Inlet Model shown 
in Figure 3. (unpublished data) 
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Figure 25. Comparisons of peak distortion factors predicted by the present 
analysis and the Melick method based on the total pressure rtns 
measurements for the subsonic inlet data set 
(unpublished data) 
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Figure 27. Comparisons of the Predicted and Measured nns Level and 

Peak Distortion Factor for the Transonic Inlet Model shown 
in Figure 4. (unpublished data) 
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Figure 28. Comparisons of peak distortion factors predicted by the present 
analysis and the Melick method based on the total pressure rms 
measurements for the transonic inlet data set 
(unpublished data) 
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Figure 29 . Continued 
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Figure 31 . Accuracy of the Present Method in Predicting rms Levels 
and Dynknic Peak Distortion Factors 
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